Fingerprints of spin-orbital physics in cubic Mott insulators: 
Magnetic exchange interactions and optical spectral weights 
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The temperature dependence and anisotropy of optical spectral weights associated with different 
multiplet transitions is determined by the spin and orbital correlations. To provide a systematic basis 
to exploit this close relationship between magnetism and optical spectra, we present and analyze 
the spin-orbital superexchange models for a series of representative orbital-degenerate transition 
metal oxides with different multiplet structure. For each case we derive the magnetic exchange 
constants, which determine the spin wave dispersions, as well as the partial optical sum rules. 
The magnetic and optical properties of early transition metal oxides with degenerate tig orbitals 
(titanates and vanadates with perovskite structure) are shown to depend only on two parameters, 
viz. the superexchange energy J and the ratio n of Hund's exchange to the intraorbital Coulomb 
interaction, and on the actual orbital state. In e g systems important corrections follow from charge 
transfer excitations, and we show that KCUF3 can be classified as a charge transfer insulator, while 
LaMnOs is a Mott insulator with moderate charge transfer contributions. In some cases orbital 
fluctuations are quenched and decoupling of spin and orbital degrees of freedom with static orbital 
order gives satisfactory results for the optical weights. On the example of cubic vanadates we 
describe a case where the full quantum spin-orbital physics must be considered. Thus information 
on optical excitations, their energies, temperature dependence and anisotropy, combined with the 
results of magnetic neutron scattering experiments, provides an important consistency test of the 
spin-orbital models, and indicates whether orbital and/or spin fluctuations are important in a given 
compound. [Published in: Phys. Rev. B 72, 214431 (2005).] 

PACS numbers: 75.30.Et, 78.20.-e, 71.27. +a, 75.10.-b 



I. SUPEREXCHANGE AND OPTICAL 
EXCITATIONS AT ORBITAL DEGENERACY 



The physical properties of Mott (or charge transfer) 
insulators are dominated by large on-site Coulomb inter- 
actions oc U which suppress charge fluctuations. Quite 
generally, the Coulomb interactions lead then to strong 
electron correlations which frequently involve orbitally 
degenerate states, such as 3c? (or Ad) states in transition 
metal compounds, and are responsible for quite com- 
plex behavior with often puzzling transport and mag- 
netic properties^ The theoretical understanding of this 
class of compounds, with the colossal magnetoresistance 
(CMR) manganites as a prominent example)^ has sub- 
stantially advanced over the last decade^ after it became 
clear that orbital degrees of freedom play a crucial role in 
these materials and have to be treated on equal footing 
with the electron spins, which has led to a rapidly devel- 
oping field — orbital physics £ Due to the strong onsite 
Coulomb repulsion, charge fluctuations in the undoped 
parent compounds are almost entirely suppressed, and 
an adequate description of these strongly correlated in- 
sulators appears possible in terms of superexchange^ At 



orbital degeneracy the superexchange interactions have a 
rather rich structure, represented by the so-called spin- 
orbital models, discovered three decades agOji& and ex- 
tensively studied in recent years.V'VMVVVVW 8 

Although this field is already quite mature, and the 
first textbooks have already appeared^*^ it has been 
realized only recently that the magnetic and the opti- 
cal properties of such correlated insulators with partly 
filled d orbitals are intimately related to each other, be- 
ing just different experimental manifestations of the same 
underlying spin-orbital physicsi^fliSi While it is clear that 
the low- energy effective superexchange Hamiltonian de- 
cides about the magnetic interactions, it is not imme- 
diately obvious that the high- energy optical excitations 
and their partial sum rules have the same roots and may 
be described by the superexchange as well. In fact, this 
interrelation between the magnetic and the optical prop- 
erties makes it necessary to reanalyze the spin-orbital su- 
perexchange models, and to extract from them important 
constraints imposed by the theory on the system param- 
eters. We will show that also the opposite holds — some 
general rules apply for the magnetic interactions in the 
correlated insulators with degenerate (or almost degener- 
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ate) orbitals, and therefore the magnetic measurements 
impose constraints on any realistic theory. At the same 
time, we shall argue that such experiments provide very 
useful information concerning the orbital order (00) and 
the strength of quantum fluctuations in a given com- 
pound, which can next be employed to interpret other 
experiments, including the optical spectroscopy. 

The phenomena discussed in the present paper go well 
beyond the more familiar situation of a Mott insulator 
without orbital degeneracy, or when the orbital degen- 
eracy is lifted by strong Jahn- Teller (JT) distortions as 
for example in the high T c cuprate superconductors. In a 
Mott insulator the optical conductivity is purely incoher- 
ent, and the optical response is found at energies which 
exceed the optical gap. When orbital degrees of freedom 
are absent, the optical gap is determined by the intraor- 
bital Coulomb interaction element U. Naively, one might 
expect that the high-energy charge excitations at energy 
~ U , which contribute to the optical intensities, are unre- 
lated to the low-energy magnetic phenomena. However, 
both energy scales are intimately related as the superex- 
change follows from the same charge excitations which 
arc detected by the optical spectroscopy. The prominent 
example of this behavior is the nondegenerate Hubbard 
model, where the virtual high-energy excitations deter- 
mine the superexchange 6 * energy J — it decides, together 
with spin correlations, about the spectral weight of the 
upper Hubbard band at half-fillingiSSiS 3 . When tempera- 
ture increases to an energy scale ~ J, the spin correla- 
tions are modified and the total spectral weight in the 
optical spectroscopy follows these changes^ 

The superexchange models for transition metal per- 
ovskites with partly filled degenerate orbitals have a 
more complex structure than for nondegenerate orbitals 
and allow both for antifcrromagnetic (AF) and for ferro- 
magnetic (FM) superexchangeii^ These different contri- 
butions to the superexchange result from the multiplet 
structure of excited transition metal ions which depends 
on the Hund's exchange Jh and generates a competi- 
tion between high-spin and low-spin excitations. The ex- 
change interactions arc then intrinsically frustrated even 
on a cubic lattice, which enhances quantum effects both 
for e g jiifiiii and for ti g systemsii^i^ This frustration is 
partly removed in anisotropic AF phases, which break 
the cubic symmetry and effectively may lead to dimen- 
sionality changes, such as in j4-type AF phase realized in 
LaMn03, or in C-type AF phase in LaV03. 

While rather advanced many-body treatment of the 
quantum physics characteristic for spin-orbital models 
is required in general, we want to present here certain 
simple principles which help to understand the heart of 
the problem and give simple guidelines for interpreting 
experiments and finding relevant physical parameters of 
the spin-orbital models of undoped cubic insulators. We 
will argue that such an approach based upon classical 00 
is well justified in many known cases, as quantum phe- 
nomena are often quenched by the Jahn- Teller (JT) cou- 
pling between orbitals and the lattice distortions, which 



are present below structural phase transitions and induce 
00 both in spin-disordered and in spin-ordered phases. 25 
However, we will also discuss the prominent example of 
LaV03, where assuming perfect 00 or attempts to de- 
couple spin and orbital fluctuations, 26 fail in a spectac- 
ular way and give no more than a qualitative insight 
into certain limiting cases. Significant corrections due 
to quantum phenomena that go beyond such simplified 
approaches are then necessary for a more quantitative 
understanding. 

In the correlated insulators with partly occupied de- 
generate orbitals not only the structure of the superex- 
change is complex, but also the optical spectra ex- 
hibit strong anisotropy and temperature dependence 
near the magnetic transitions, as found in LaMn03j2L2£ 
the cubic vanadates LaV03 and YVQ 3 ) 29 i 30 and in the 
ruthenatcs. 31 In all these systems several excitations con- 
tribute to the excitation spectra, so one may ask how 
the spectral weight redistributes between individual sub- 
bands originating from these excitations. The spectral 
weight distribution is in general anisotropic already when 
00 sets in and breaks the cubic symmetry, but even more 
so when A-type or C-type AF spin order occurs below the 
Neel temperature. 

The effective spin-orbital models of transition metal 
oxides with partly filled degenerate orbitals depend in a 
characteristic way upon those aspects of the electronic 
structure which decide whether a given strongly corre- 
lated system can be classified as a Mott insulator or as a 
charge transfer (CT) insulator. As suggested in the orig- 
inal classification of Zaanen, Sawatzky and Allen^ the 
energy of the d — p CT excitation A has to be compared 
with the Coulomb interaction U — if U < A, the first 
excitation is at a transition metal ion and the system is a 
Mott insulator, otherwise it is a CT insulator. Both are 
strongly correlated insulators, yet in one limit the domi- 
nant virtual excitations are of d — d type, whereas in the 
other limit they are of p — d type. One may consider 
this issue more precisely by analyzing the full multiplet 
structure, and comparing the lowest excitation energy 
(to a high-spin configuration) at a transition metal ion, 
£ns = U — 3Jh, with that of the lowest CT excitation 
(of energy A) between a transition metal ion and a ligand 
ionm^ Thus we argue that one can regard a given per- 
ovskite as a charge transfer insulator if ens > A, and as 
a Mott-Hubbard insulator if ens < A. By analyzing these 
parameters it has been suggested that the late transition 
metal oxides may be classified as CT insulators i In this 
case important new contributions to the superexchange 
ar i sej 34 l 35 i 3£ ca j]_ed below A (charge transfer) terms, as we 
shall discuss for two e g systems: KCUF3 and LaMn03. 

A central aim of this paper is to provide relatively sim- 
ple expressions for the magnetic exchange constants and 
for the optical spectral weights that can be used by exper- 
imentalists to analyze and compare their spin wave data 
with optical data. While the full spin-orbital models are 
rather complex, they are nevertheless controlled by only 
very few physical parameters: (i) the superexchange con- 
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stant J, (ii) the normalized Hund's exchange 77, and (iii) 
the charge transfer parameter R. There are two distinct 
ways to determine these effective parameters: either (i) 
from the original multiband Hubbard model, or (ii) from 
experimental spin wave and/or optical data. Here the 
second approach is of particular interest because the si- 
multaneous analysis of magnetism and optics provides a 
subtle test of the underlying model. 

The paper is organized as follows. In Sec. II we in- 
troduce the generic structure of the low-energy effective 
Hamiltonian in a correlated insulator with orbital de- 
generacy, and discuss its connection with the optical ex- 
citations at high-energy. This general formulation pro- 
vides the important subdivision of a given spin-orbital 
model which is necessary to obtain the partial spectral 
weights for individual multiplet transitions. In the re- 
maining part of the paper we concentrate on some se- 
lected cubic perovskites and demonstrate that this gen- 
eral formulation allows one to arrive at a consistent in- 
terpretation of the magnetic and optical experiments in 
these correlated insulators using the superexchange in- 
teractions (Sees. Ill- VI), and to deduce the parameters 
relevant for the theoretical model from the experimental 
data, whereever available. We start in Sec. Ill with the 
simplest spin-orbital model for e g holes in KCUF3, and 
demonstrate that this system is in the CT regime which 
changes the commonly used picture of superexchange in 
this system in a qualitative way. Next we present and 
analyze the spin-orbital model with e g orbital degrees 
of freedom for the undoped manganite LaMnOs in Sec. 
IV. Here we show that in this case much smaller con- 
tributions arise from the CT processes, and the system 
is already in the Mott-Hubbard regime of parameters, 
which explains the earlier success of a simplified effective 
model based entirely on d — d excitations and sufficient 
for a semiquantitative understanding. This justifies our 
approach to the early transition metal perovskites with 
t2 g degrees of freedom, titanates in Sec. V and vanadates 
in Sec. VI, which we treat as Mott-Hubbard insulators. 
For all these systems we analyze the magnetic exchange 
interactions and the optical spectral weights, depending 
on the nature of the spin correlations in the ground state. 
The paper is concluded in Sec. VII, where we provide a 
coherent view on the magnetic and the optical phenom- 
ena and summarize the experimental constraints on the 
model parameters. 

II. GENERAL FORMALISM 

We consider here effective models with hopping ele- 
ments between transition metal ions, 

Here Ei a are orbital energies, and tiajp are effective hop- 
ping elements via ligand orbitals — they depend on the 
type of considered orbitals as discussed in Refs. |37] and 



l38l The energy scale for the hopping is set by the largest 
hopping element t: the {dda) element in case of e g sys- 
tems, and the (ddn) element when only it bonds are con- 
sidered in systems with degenerate and partly filled ti g 
orbitals. For noninteracting electrons the Hamiltonian 
Hq would lead to tight-binding bands, but in a Mott insu- 
lator the large Coulomb interaction U suppresses charge 
excitations in the regime of U ^> t, and the hopping ele- 
ments can only contribute via virtual excitations, leading 
to the superexchange. 

The superexchange in the 3d cubic systems with orbital 
degeneracy is described by spin-orbital models, where 
both degrees of freedom are coupled and the orbital state 
(ordered or fluctuating) determines the spin structure 
and excitations, and vice versa. The numerical and an- 
alytical structure of these models represents a fascinat- 
ing challenge in the theory, as it is much more complex 
than that of pure spin models. The spin-orbital models 
have been derived before in several cases, and we refer 
for these derivations to the original literature>ii*iiiSii& 
They describe in low energy regime the consequences of 
virtual charge excitations between two neighboring tran- 
sition metal ions, # c£™ +1 <i™ -1 , which involve an 
increase of energy due to the Coulomb interactions. Such 
transitions are mediated by the ligand orbitals between 
the two ions and have the same roots as the superex- 
change in a Mott insulator with nondegenerate orbitals 6 
at U ^ t — thus the resulting superexchange interac- 
tions will be called U terms. The essential difference 
which makes it necessary to analyze the excitation ener- 
gies in each case separately is caused by the existence of 
several different excitations. Their energies have to be 
determined first by analyzing the eigenstates of the local 
Coulomb interactions, 

Hint = U^^n ia ^n ia i + (U a (3 ~ ^Jaf^jniaTUp 

ia. i,a</3 

i,a<f3 

— 2 JafiSia ■ S^, (2.2) 

i,a<f3 

with a = —a, which in the general case depend on the 
three Racah parameters A, B and which may be 
derived from somewhat screened atomic values. While 
the intraorbital Coulomb clement 

U = A + AB + 3C, (2.3) 

is identical for all 3d orbitals, the interorbital Coulomb 
and exchange elements, U a p and J a p, are anisotropic and 
depend on the involved pair of orbitals; the values of 
J a j3 are given in Table [J The Coulomb and exchange 
elements are related to the intraorbital element U by a 
relation which guarantees the invariance of interactions 
in the orbital space, 

U = U al3 + 2J af3 . (2.4) 
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In cases where only the orbitals of one type (e g or t2 g ) 
are partly filled, however, as e.g. in the titanates, vana- 
dates, or copper fluorides, all relevant exchange elements 
J a p are the same (see Table HJ and one may use a sim- 
plified form of onsite interactions^ 



H 



(0) 



U 



/ a 

Jh V U in A 



i,Q</3 



/ , y^ial u ial u iPl u i/3| ^ "i^f ««i 
i,Q</3 



2J^ ^ Sj Q ■ S^. 

i,a</3 



(2.5) 



with only two parameters: the Coulomb element U (|2.3|) 
and a Hund's exchange element Jh, being AB + C for 
e g and 3_B + C for systems, respectively. We em- 
phasize that in the general case when both types of or- 
bitals are partly filled (as in the manganites) and both 
thus participate in charge excitations, the Hamiltonian 
(|2.5() is only approximate, and the full excitation spec- 
tra of the transition metal ions^ which follow from Eq. 
(|2.2|) have to be used instead. A few examples of spectra 
for d^d" 1 # dT^dT -1 charge excitations at transition 
metal ions are shown in Fig. [T] As a universal feature, 
the high-spin excitation is found at energy U — 3Jh in 
all cases, provided that Jh is understood as Hund's ex- 
change for that partly filled manifold (e g or t2 g ) of degen- 
erate d orbitals which participate in charge excitations. 
The structure of the excited states depends on the partly 
occupied orbitals^i and on the actual valence m — the 
distance between the high-spin and low-spin excitations 
increases with the number of electrons for m < 5 (holes 
for m > 5). 

At orbital degeneracy the superexchange which con- 
nects ions at sites i and j along the bond (ij) involves 
orbital operators which depend on the bond direction. 
Therefore, we introduce the index 7 = a, b, c to label the 
three a priori equivalent directions in a cubic crystal. 
In order to analyze the consequences of each individual 
charge excitation n that contributes to the superexchange 
in a given transition metal compound with degenerate d 
orbitals, we shall use below a general way of writing the 
effective low-energy Hamiltonian as a superposition of 



TABLE I: On-site interorbital exchange elements J a/ 3 for 3d 
orbitals as functions of the Racah parameters B and C (for 
more details see Ref. KM). 
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FIG. 1: Energies of dfdf -> charge excitations 

in selected cubic transition metal oxides, as obtained from 
Eq. for: (a) e g excitations of Cu 3+ (d 8 ) and Mn 2+ (d 5 ) 

ions [in the d 5 case the spectrum— was obtained from Eq. 
lEU]; (b) t 2g excitations of Ti 2+ (d 2 ) and V 2+ (d 3 ) ions. The 
splittings between different states are due to Hund's exchange 
element Jh which refers to a pair of e g electrons in (a), and 
to a pair of t 2g electrons in (b), respectively. 



such individual terms on each bond (ij), 

" (<3>||7 



(2.6) 



with the energy unit [absorbed in individual H^\ij) 
terms] given by the superexchange constant, 



J 



4*2 
~U~' 



(2.7) 



It follows from d — d charge excitations with an effective 
hopping element t between transition metal ions, and is 
the same as that obtained in a Mott insulator with non- 
degenerate orbitals in the regime of U 3> t& While U 
is the uniquely defined on-site intraorbital Coulomb ele- 
ment l|2.H|l . increasing upon going from Ti to Cu along 
the transition metal series^S^ the definition of the hop- 
ping t between two nearest neighbor transition metal ions 
depends on the system^ If degenerate e g orbitals are in- 
volved, it is the effective (dda) hopping element for a er- 
bond which involves p a orbitals on the intervening ligand 
ion (e.g. for the hopping between two directional 3z 2 — r 2 
states along the c axis), while for the systems with degen- 
erate t2 g orbitals it stands for the effective (ddir) hopping 
element due to tt bonds which involve p^ orbitals on the 
ligand ion. 

In the superexchange Hamiltonian Eq. H2.fi|l the con- 
tributions which originate from all possible virtual ex- 
citations d™d™ ^ d- n+1 d™ _1 just add up to the total 
superexchange interaction, in which the individual terms 
cannot be distinguished. Yet each of these excitations 
involves a different state in the multiplct structure of at 
least one of the transition metal ions, i.e., either in the 
d m+ or in the d™" 1 configuration or in both, depending 
on the actual process and on the value of m. As pointed 
out elsewhere^ the same charge excitations contribute 
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to the optical conductivity, and here they appear at dis- 
tinct energies, thus revealing the multiplet structure of 
the excited transition metal ions. Moreover, they convey 
a rich and characteristic temperature dependence to the 
optical spectrum, determined by the temperature varia- 
tion of the spin-spin and orbital-orbital correlations. We 
emphasize that it is therefore important to analyze the 
various multiplet excitations separately, as they depend 
on these correlations in a different way, and will also con- 
tribute to a quite different temperature dependence, as we 
show in this paper on several examples. 

As we will see in more detail below, the generic struc- 
ture of each such individual contribution is for a bond 
(ij) given by 

H<?\ij) = (a n + b n S l -S J )Q^(r l: r j ) 

+b n Q^(r i ,r j )S i -§ j , (2.8) 

where the orbital dependence of the superexchange is de- 
scribed by means of orbital projection operators Qn^ 
which are expressed in terms of components of orbital 
pseudospin T = 1/2 operators at sites i and j. The 
coefficients a n and 6 n , which measure the strength of 
the purely orbital part and of the spin-and-orbital part 
of the superexchange, respectively, follow from second- 
order perturbation theory involving the charge excitation 
ii. In the present case of perovskites, where the bond be- 
tween two transition metal ions through the ligand ion 
(F or O) connecting them is close to linear (180°), the 
coefficients a n and b n are of similar magnitude (in con- 
trast to the situation in layered compounds like LiNi02 
with 90° bonds where the purely orbital interaction is 
stronger by an order of magnitude than the spin-and- 
orbital interaction^) . 

Here we consider systems having cubic symmetry at 
high temperature. Yet at low temperature this symme- 
try is frequently spontaneously broken — usually driven 
by the joint effect of (i) the orbital part of the superex- 
change interaction, and (ii) the JT coupling of the same 
degenerate (and therefore JT active) 3d orbitals to lattice 
modes. The result is the simultaneous onset of a macro- 
scopic lattice distortion and of 00, i.e., a cooperative JT 
effect. At temperatures well below the transition temper- 
ature T s of this combined structural and orbital phase 
transition, the 00 is effectively frozen. The remaining 
superexchange interactions between the spins may then 
be obtained by replacing the orbital projection operators 
in Eq. 12.611 by their expectation values, 

Qn\n,Tj) — (QtffaZj)) = Q^m, (t,-». (2.9) 

Obviously, this leads to anisotropic magnetic interac- 
tions, 

ff, = j^^MQl 7) )Si'4 (2.io) 

n <y)ll7 



which will in general induce a further magnetic phase 
transition at lower temperature. It is noteworthy that 
in this situation the spin degrees of freedom get decou- 
pled from the orbital degrees of freedom, although the 
purely orbital (a n ) and spin-and-orbital (b n ) superex- 
change terms are of similar strength. Responsible for 
this behavior is the JT contribution to the structural 
phase transition, which enhances T s above the value it 
would have if the transition were driven by orbital su- 
perexchange alone. 

Starting from the microscopic spin-orbital superex- 
change models, we will analyze the effective spin mod- 
els which arise after such a symmetry breaking at low 
temperature. Rewritten from Eq. (|2.1U|) . they are of the 
generic form 

h s = j c ^ ■ Sj + j ab Si ■ 4 ( 2 - n ) 

{ij)c {ij)ab 

with two different effective magnetic exchange interac- 
tions: J c along the c axis, and J ab within the ab planes. 
The latter J ab interactions could in principle still take 
two different values in case of incquivalent lattice dis- 
tortions (caused, e.g., by octahedra tilting or pressure 
effects) making the a and b axes inequivalent, but we in- 
tend to limit the present analysis to idealized structures 
with these two axes being equivalent. We shall show that 
the spin-spin correlations along the c axis and within the 
ab planes, 

Sc = (Si ■ Sj)c, S ab = (Si ■ Sj)ab, (2-12) 

next to the orbital correlations, play an important role 
in the intensity distribution in optical spectroscopy. 

The spectral weight in the optical spectroscopy is de- 
termined by the kinetic energy^ and reflects the onset 
of magnetic order^ii and/or orbital order As shown 
by Ahn and Millis^ in the weak coupling regime one 
can analyze the total spectral weight in optical absorp- 
tion using the Hartree-Fock approximation for the rele- 
vant tight-binding Hamiltonian. In a correlated insula- 
tor the electrons are almost localized and the only ki- 
netic energy which is left is associated with the same 
virtual charge excitations that contribute also to the su- 
perexchange. Therefore, we will discuss here the indi- 
vidual kinetic energy terms K^C\ which can be deter- 
mined from the superexchange (|2.6|) using the Hellman- 
Feynman theorem^ 

K<?~>=-2(H<?\ij)). (2.13) 

For convenience, we define the as positive quanti- 

ties. Making use of the generic form of the superexchange 
contribution H^\ij) given by Eq. i|2.8[l . and assuming 
as above that spin and orbital degrees of freedom are de- 
coupled in the temperature range of interest, we obtain 

K™ = -2j(a n + b n (St ■ S^) (Q^(n,T 3 )). (2.14) 
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Each term i|2.13|l originates from a given charge 

excitation n for a bond (ij) along axis 7. These terms 
are directly related to the partial optical sum rule for 
individual Hubbard bands, which reads 20 



aoh 



r (7) 



(u))du> 



2 " ' 



(2.15) 



where ali\u>) is the contribution of band n to the opti- 
cal conductivity for polarization along the 7 axis, ao is 
the distance between transition metal ions, and the tight- 
binding model with nearest neighbor hopping is implied. 
Comparison with Eq. Ij2.14|) shows that the intensity of 
each band is indeed determined by the underlying 00 to- 
gether with the spin-spin correlation along the direction 
corresponding to the polarization. 

One has to distinguish the above partial sum rule 
(|2.15(l from the sum rule for the total spectral weight 
in the optical spectroscopy for polarization along a cubic 
direction 7, involving 



(2.16) 



which stands for the total intensity in the optical excita- 
tions (due to d — d excitations) . This quantity is of less 
interest here as it has a much weaker temperature de- 
pendence and does not allow for a direct insight into the 
nature of the electronic structure. In addition, it might 
be also more difficult to resolve from experiment. 

When the low-energy excitations are of CT type, two 
holes could also be created within a 2p orbital on a ligand 
(oxygen or fluorine) ion in between two transition metal 
ions, described by d™p 6 d™ ^ d-" +1 p 4 d™ +1 processes - 
these CT contributions lead to additional superexchange 
contributions, called below A terms. While the latter 
terms can be safely neglected in Mott-Hubbard systems, 
they substantially modify the superexchange in CT in- 
sulators, and may even represent there the dominating 
contribution^^ Below we will analyze them in two sit- 
uations which involve e g degrees of freedom, viz. in the 
cubic copper fluoride KCuF 3 (Sec. Ill), and in the cubic 
manganite LaMn03 (Sec. IV), and we will show that in 
KCUF3 they represent an essential part of the superex- 
change. 



III. COPPER FLUORIDE PEROVSKITE: KCuF :i 

A. Superexchange Hamiltonian 

The simplest spin-orbital models are obtained when 
transition metal ions are occupied by either one electron 
(m = 1), or by nine electrons (m = 9); in these cases the 
Coulomb interactions l|2.5[) contribute only in the excited 
state (in the d 2 or the d 8 configuration) after a charge 
excitation d^d™ ^ between two neighboring 

ions. Here we start with the case of a single hole in the 



d shell, as realized for the Cu 2+ ions in KCUF3 with the 
d 9 configuration (m = 9). Due to the splitting of the 3d 
states in the octahedral field within the CuFg octahedra, 
the hole at each magnetic Cu 2+ ion occupies one of the 
e g orbitals. The superexchange coupling l|2.6|) is usually 
analyzed in terms of e g holes in this caseji and this has 

become a textbook example of spin-orbital physics by 

4 19 
nowi»ia 

Orbital order occurs in KCUF3 below the structural 
transition at T s — 800 K. At T < T s the structure is 
tetragonal, with longer Cu— Cu distances within the ab 
planes (d a b — 8.28 A) than along the c axis (d c = 7.85 
A) 4^ which favors strong AF interactions along the c 
axis. Below the magnetic transition at Tjv — 38 K, long- 
range magnetic order of A-type sets inj£2*£i and the or- 
dered moment is /io = 0.48/iB^ 

The superexchange between the Cu 2+ ions in KCUF3, 



H(d 9 ) =Hu(d 9 )+H A {d 9 ) 



(3.1) 



consists of two terms: the U term Hu (|2.6|) . and the CT 
term Ha- First we introduce the U term Tijj(d 9 ) follow- 
ing the general approach of Sec. [H] It originates from 
three different excitations, leading to an intermediate d s 
configuration at an excited Cu 3+ ion. Using the model 
Hamiltonian l|2.5l) to describe the Coulomb interactions 
between the e g electrons, one finds an equidistant excita- 
tion spectrum of 3 A2, 1 E ^Eg and 1 E e ) and l Ax states, 
with energies: 11 i 39 U — 3Jh, U — Jh and U + Jh, as 
shown in Fig. This excitation spectrum is exact, 

and the element Jh for a pair of e g electrons is given by 
the Racah parameters B and C (see Table I): 



J H = 45 + C. 



(3.2) 



This definition of Jh will be used for two systems with e g 
orbital degrees of freedom: for the copper fluoride KCuF 3 
(considered here), and for the manganite LaMn03 (in 
Sec. llVll . 

In what follows, we will parametrize the multiplet 
structure of the different transition metal ions by the 
ratio of the Hund's element Jh and the intraorbital 
Coulomb element U, 



>1 



Jh 
U ' 



(3.3) 



Using Eqs. I|2.6|) and (|3.3|) . one finds for each bond (ij) 
along a 7 axis (7 = a, 6, c) four contributions i— 



H 



(7) 



, (3.4) 
(3.5) 
(7) \ (3.6) 

with coupled spin and orbital operators. The coefficients, 



r(7) 



m<> = 



-)C- 

A) V4 
l\/f 
4A2 



.(7) 



1 



1 



1 



'"1 



1 — 377' 



r2 = r 3 



1-7? 



, r 4 



1 + 77 



(3.8) 
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follow from the above multiplet structure of d 8 ionsS 

As explained below, it is straightforward to under- 
stand the generic structure of the superexchange term 
Hu, given by Eqs. $T ty -$F fy . Here Si is a spin S = 1/2 
operator, and 



P 1 (ij) = S i -S j + I , 



P (ij) = ~-S v S j , (3.9) 



are the spin projection operators on triplet (S = 1) 
and singlet (S — 0) spin states for a bond (ij), so one 
recognizes the high-spin term lj3.4|l and three low-spin 
terms H3.5|) - (|3.7|l . respectively. The spin operators in 
Eqs. i|3.4fl - i|3.7fl are accompanied by orbital pseudospin 
operators oc t^ 7 ', which select the type of orbitals occu- 
pied by holes at sites i and j, and simultaneously dictate 
the allowed excited states. 

The orbital operators depend on the direction of 
a considered bond (ij), and are given by 



(3.10) 



where trf and o~f are Pauli matrices acting on the orbital 

pseudospins and the signs ± in T^ ab ^ correspond to a and 

b axis, respectively. With the help of one defines 
next the projection operators in the orbital subspace, 



Q 



.(7) 



i 



(3.11) 



For a given cubic axis 7 they project (at site i) either on 
the planar orbital |£) in the plane perpendicular to the 
7 axis, or on the orthogonal directional orbital |C) along 
this axis. For instance, in the case where 7 is the c axis, 
they project on the x 2 — y 2 orbital in the ab plane, and 
on the directional 3z 2 — r 2 orbital along the c axis. 

Using the projection operators (|3.11() . the orbital de- 
pendence in Eqs. I|3.4|) - (|3.7|) becomes transparent. First 
of all, (| — r^ 7 ' 'rj 7 ' ') in Eq. 1|3.4|) accompanies the high- 
spin 3 A2 excitation as this state may occur only when a 
pair of orthogonal orbitals is occupied at sites i and j, 
described formally by a superposition of two possibilities, 
\{QiC,Qjt, + Qi^Qjc)- In contrast, the operator Qi^QjQ 
selects two orbitals oriented along the bond (ij) for the 
high-energy low-spin l A\ excitation, see Eq. (|3.7|) . Fi- 



nally, the second Hy' (E31> and third (E2IJ) term 

correspond to the doubly degenerate low-spin 1 E state 
which consists of two singlet excitations: (i) an interor- 
bital singlet with two different orbitals occupied ( 1 i? e ), 
and (ii) a double occupancy within a directional orbital at 
either site ( 1 £'e) — these two excitations have thus quite 
different orbital dependences, identical with those of the 
3 A2 and the ^1 excitation, respectively. The sum over 
all the terms Hn \ with n = 1, • • • , 4, gives the simplest 
version of the spin-orbital model for the cubic copper flu- 
oride KCUF3 with degenerate e g orbitals. Its derivation 
and more details on the classical phase diagram can be 
found in Ref. Illl 



By considering further the electronic structure of 
KCuF 3 , one can elucidate the role played by the CT part 
Ha in the superexchange Hamiltonian 13.1fl . By analogy 
with the CuC>2 planes of the high-temperature supercon- 
ductors, where the CT processes give the dominating 
contribution to the AF superexchange interaction's^ 
one expects that they are also important for a cubic cop- 
per (II) fluoride and modify the superexchange in KCUF3. 
The CT term, 



H A (d 9 



JR 



E 

(*j)ll7 



[Si- 



Si 



J7) 



with the coefficient 



R = 



2U 



2A + U V 



„(7) 



(3.12) 



(3.13) 



resulting from the two-hole charge excitation at a com- 
mon neighboring 2p a orbital of a fluorine ion in between 
two copper ions, in the process d?p^d® dj°p%^d^ . 
As a double hole excitation is generated at a single bond- 
ing orbital 2p a within each Cu— F— Cu unit, this term 
is necessarily AF. Two holes can move to fluorine from 
two neighboring Cu ions only if both of them occupy a 
directional e g orbital \Q, oriented along the considered 
bond (e.g., 3ir 2 — r 2 orbitals along the a axis), being the 
simplest CT term discussed recently by Mostovoy and 
KhomskiiiSi Therefore this process leads to the same or- 
bital dependence in Eq. I|3.12|l as the low-spin 1 Eg and 
M.1 excitations which involve double occupancies of di- 
rectional |C) orbitals. 



B. Spin exchange constants and optical intensities 

A characteristic feature of spin-orbital superexchange 
models with e g orbital degrees of freedom is the presence 
of the purely orbital interactions in Eqs. (|3.4|I - H3.7|I and 
(|3.12|l . which favor particular type of occupied orbitals. 
LDA+U calculations^^ have indicated that such purely 
electronic interactions would already drive the instability 
towards the C-type 00 (C-OO) phase, with alternating 
orbitals in the ab planes, and repeated orbitals along the 
c axis, which induces FM spin exchange in the ab planes, 
and strong AF exchange between the planes. Experimen- 
tally, this 00 sets in below the structural transition at 
T s ~ 800 KfSi i.e., at much higher temperature than the 
characteristic energy scale of the magnetic excitations, 57 
suggesting that the JT effect plays an important role in 
this instability. This observation is consistent with the 
large difference between T s and the Neel temperature 
Tjv ~ 38 Kj52 the latter being controlled by the magnetic 
part of the superexchange, and thus the orbital correla- 
tions decouple from the spin-spin correlations. This mo- 
tivates one to analyze the dependence of the magnetic 
exchange interactions and of the optical spectral weights 
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FIG. 2: Schematic representation of e g orbitals as obtained 
for different angle 6 in Eq. 13. 1511 . Pairs of orthogonal orbitals, 
forming a basis in the orbital space, differ by angle 9 — n. 
Directional e 9 orbitals \() along a and 6 axes are obtained from 
the {|3z 2 — r 2 }, \x 2 — y 2 )} basis 113. 1611 by the transformation 
lTH5t with angle 6 = ±2tt/3. Dashed lines for 6 = ±tt/3 
indicate the possible OO of \y 2 — z 2 )/\x 2 — z 2 ) orbitals occupied 
by holes in KCUF3 , suggested for T < Tn in Ref . 0. 



on the type of OO stabilized below the structural transi- 
tion. 

Here we are interested in the low temperature range of 
T < 500 K, so we assume perfect OO given by a classical 
ansatz for the ground state, 

i*o> = n i^>* n I**)* ( 3 - w ) 

with the orbital states, \0a)i and |0b)j, characterized by 
opposite angles {9 a — —9 b) and alternating between two 
sublattices A and B in the ab planes. The orbital state 
at site i: 

\e)=co S Q\z)+smQ\x), (3.15) 

is here parametrized by an angle 9 which defines the am- 
plitudes of the orbital states 

\z) = (3z 2 - r 2 )/V6, \x) = (x 2 - y 2 )/V2, (3.16) 

being a local e g orbital basis at each site. This and other 
equivalent orbital bases are shown schematically by pairs 
of solid and dashed lines (corresponding to pairs of or- 
bitals {|£), \9 + 71-)}) in Fig. H The OO state specified in 
Eq. (pTlljl is thus defined by: 

\6 A )i = cos(^-J|z) J +sin(^-J|a;) i , 
\0b) 3 = cosQl^ -sinQlx),, (3.17) 
with 9 A = 9 and 9 B = -9. 



The magnetic superexchange constants J a b and J c in 
the effective spin model (|2.11ll are obtained by decou- 
pling spin and orbital variables and next averaging the 
orbital operators in the spin-orbital model (|3.1|l over the 
classical state |<&o) as given by Eq. (|3.14l) . The relevant 
averages are given in Table ^ and they lead to com- 
pact expressions for the superexchange constants in Eq. 

J c = -jj - ri sin 2 9 +(r 2 +r 3 )(l + cos 9) 

+ (r 4 + 2i?)(l + cos6») 2 }, (3.18) 

Jab = ijj-rr^+sin 2 ^ +(r 2 + r 3 ) (l-i cos 6») 

+ (r 4 + 2i?)(i-cos6l) 2 }, (3.19) 

which depend on three parameters, viz. J l|2.7|l . 7/ l|3.3|) 
and R l|3.13|l . and on the OO (|3.17|) specified by the or- 
bital angle 9. It is clear that the FM term cx r\ competes 
with all the other AF low-spin terms. Nevertheless, in 
the ab planes, where the occupied e g orbitals alternate, 
the large FM contribution (when sin 2 9 ~ 1) still makes 
the magnetic superexchange J a b weakly FM (J a b ^ 0), 
while the stronger AF superexchange along the c axis 
{J c ^ I J a b\) favors quasi one-dimensional (ID) spin fluc- 
tuations. 

By considering the superexchange model, one can de- 
rive as well the pure orbital interactions which stabilize 
the OO. The superexchange interactions are anisotropic 
below the structural transition at T s . In contrast, at 
sufficiently high temperature T > T s , when also spin 
correlations may be neglected, one finds isotropic orbital 
interactions, 

Jc = J°al h = ^(3ri -u- 2R) , (3.20) 

which multiply t\ Tj for each bond, contributing to an 
orbital instability towards alternating G-type OO, while 
actually C-type OO is observed below T s . It is thus 
clear from experiment that this instability cannot be of 
purely electronic origin, and that, similarly to what is 
the case in LaMn03j£& it is supported by the lattice. In 
fact, although it has been argued that the OO is caused 



TABLE II: Averages of the orbital projection operators 
standing in the spin-orbital interactions in Eqs. I3.4H - 13.7II 
for the C-type (or G-type) OO of occupied e g orbitals which 
alternate in ab planes, as given by Eqs. 13.171 . Nonequivalent 
cubic directions are labelled by 7 = ab, c. 



7 ab c 

(i_ r W r W) i(f+sin 2 0) isin 2 

(G+^XH^)) KHcose) 2 i(l-cose) 2 
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primarily by the superexchange^ the electronic interac- 
tions 1)3. 20[) predict that T s ~ 0.1J, and not T s ~ J, as 
observed. Note also that as soon as the AF spin cor- 
relations develop along the c axis, one finds anisotropic 
orbital interactions, J^ b > JJ, which amplifies the ongo- 
ing symmetry breaking in the tetragonal phase. 

The spectral weights of the optical subbands also follow 
from the superexchange processes, and are determined 
from the effective Hamiltonian (|3.1|> by the general re- 
lations given by Eqs. I|2.13[) and l|2.15[) . Following the 
excitation spectrum of Fig. ^ one finds optical absorp- 
tion at three different energies (the degeneracy of the 1 E 
state is not removed), so we label the respective kinetic 
energy contributions ifi 7 ' by n — 1,2,3. They are de- 
termined at low temperature T < 500 K by rigid C-type 
00 i|3.17fl . with the classical averages of the orbital op- 
erators given in Table [H] So one finds for polarization 
along the c axis 

k[ c) = ijn (| +s c ) sin 2 9, (3.21) 

4 C) = ^J(r 2 + r 3 )(^-s c )(l + cos 0), (3.22) 

K { 3 c) = ijr 4 (i - s c )(l + cos6») 2 , (3.23) 
and for polarization in the ab plane 

K[ ab) = ijn (! + *.*)(!+ Bin* *), (3.24) 

K ( 2 ab) = ij(r 2 +r 3 )(i- 5ah )(l-icos0), (3.25) 

*t 6) = ijr 4 (i- Safc )(i-cos0) 2 . (3.26) 

Similar to the exchange constants J a b and J c , the kinetic 
energies depend on the multiplet structure described by 
two parameters, viz. J (|2.7|) . r\ (|3.3|l . and on the 00 
(|3.17l) specified by the angle 9. Note that they depend 
on these parameters also indirectly, since the spin-spin 
correlations are governed by J a b and J c as well. We an- 
alyze this dependence in Sec. IIII Dl 

C. Magnetic interactions in KCuFa 

In order to apply the above classical theory to KCUF3 
we need to determine the microscopic parameters which 
decide about the superexchange constants, given by Eqs. 
(|3.18|l and (|3.19|l . In principle, if the optical data would 
also be available, with this experimental input one would 
be able to fix the values of the relevant parameters J and 
?7, and the orbital angle 9. Having only magnetic mea- 
surements, we give here an example of another approach 
which starts from the microscopic parameters for the lo- 
cal Coulomb interaction and Hund's element suggested 
by the electronic structure calculations performed within 
the LDA+U methodt^ U = 7.5, J H = 0.9 eV — they 
lead to T) — 0.12. Note that these parameters are some- 
what smaller than the values U = 8.96 and Jh = 1.19 




6 (degrees) 

FIG. 3: (Color online) Exchange interactions (13.1811 : J c 113.1811 
and Jab (13. 1911 for the copper fluoride model as functions of 
the orbital angle 9 which describes the OO [see Eq. I|3.14jl ]. 
The dashed lines show the U term alone, while the solid lines 
include the CT contributions as well. The OO induced by the 
JT distortions (6 ~ 70°) is indicated by dotted line. Param- 
eters: J ~ 33.3 meV, r\ = 0.12, R = 1.2. 



eV deduced for Cu 2+ ions in the Cu0 2 planes of the 
high-temperature superconductors by Grant and McMa- 
han using the fixed charge method^ but we believe that 
they reflect better the partly screened interactions within 
CuFg units. We are not aware of any estimation of the 
remaining microscopic parameters until now, but taking 
into account the expected contraction of the 2p wavefunc- 
tions by going from 2 ~ to F - ions, we argue that t p d is 
reduced, while A and U p could be similar to their respec- 
tive values for Cu0 2 planes. 5 ^ Therefore, it is reasonable 
to adopt: t pd — 1.0, A = 4.0, and U p = 4.5 eV. Note that 
although the values of t p d = 1.0, A and U p could not be 
really estimated, in the present approach they are not 
independent parameters; also only a linear combination 
of A and U p enters Eq. H3.13|) . so a change in the value 
of U p could to some extent compensate a modified value 
of A. The present parameters lead to J ~ 33.3 meV and 
R=l.2. 

Consider now the OO of the occupied orbitals (by 
holes) in KCUF3. Recent resonant x-ray scattering exper- 
iments suggest that both sublattices are equivalent, with 
9 A = -Ob = 9 in Eq. I|3.14|l . 5 - but the precise shape 
of the occupied orbitals in KCUF3 remains unresolved. 
There are different views concerning the type of orbitals 
that participate in the OO state. On the one hand, it 
is believed that the orbital angle 9 should be close to 
the angle 6*jt — 70° (~ 0.397r), as given by cos 9 ~ |, 
which follows from the local lattice distortions M> On the 
other hand, the electronic interactions in the symmetry 
broken A-AF phase below Tm would favor instead alter- 
nating (y 2 — z 2 )/(x 2 — z 2 ) orbitals, 7 with 6> SE = 60°. In 
reality, one expects rather a certain compromise between 
the electronic interactions for finite spin correlations and 
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those induced by the lattice. Thus, in the present study 
of the magnetic exchange constants and optical spectral 
weights we shall consider a range of possible values of 
60° < 9 < 90°, focusing in particular on the above val- 
ues favored by the above individual terms in the effective 
Hamiltonian. 

First, we demonstrate that the model Eq. 1)3. 1|) is ca- 
pable of describing the experimentally observed exchange 
constants J° xp ~ 35 meV and J™ p ~ -2 meVM. Re- 
markably, the value of |J°^ P | is smaller by more than 
one order of magnitude than J° xp , being some challenge 
for the theoretical model. Consider first the values of J c 
and J a b for varying angle 9 which tunes the 00 (Fig. 
El). When only the U part of the superexchange 7iu is 
considered, one finds J c ~ 15 meV and J ao ~ —2 meV. 
The FM term oc r\ gives the largest contribution for the 
ab planes, which follows from the alternation of hole or- 
bitals in the ab planes, being close to the planar orbitals 
(x 2 — z 2 )/(y 2 — z 2 ) (9 = 60°) suggested early on by Kugel 
and Khomskii, 7 but is partly compensated by the AF 
terms. While J a b is only weakly depending on 9 near 
this type of OO, J c decreases steadily with increasing 9 
(Fig. |3) , as the overlap between the orbitals occupied by 
holes along the c axis decreases when the amplitude of 
the \z) states is reduced. One finds that using only the U 
term in the superexchange, rather extreme parameters, 
such as U < 4 eV and Jh ~ 0.3 eV (with the present 
value of t) , would have to be assumed to reproduce the 
experimental values of J c and J ao . 

The CT term (jXT^I with R = 1.2 enhances the AF 
interaction J c by a factor larger than two but hardly 
changes J a b (see Fig. |3J). Only then J c comes close to the 
experimental value J° xp ~ 35 meV^ As in the U model 
(at R = 0), the value of J c decreases with increasing 9, 
and there is no serious difficulty to fit the parameters in 
order to obtain a reasonable agreement with experiment, 
once the value of the orbital angle 9 would be known. As 
an illustrative example we show the results obtained with 
the present parameters in Fig. [3] — one finds J c ~ 32 
meV for the OO which agrees with the lattice distortions 
(0j T ^ 70°), and J c ~ 40 meV for the Kugel-Khomskii 
(x 2 - z 2 )/(y 2 - z 2 ) orbitals (0 SE = 60°). These results 
demonstrate that the CT superexchange term plays an 
essential role in KCUF3 and so this system has to be 
classified as a charge transfer insulator. 

Remarkably, the value of J a b is almost unaffected by 
the CT term (Fig. |3J). This is due to the alternating OO 
in the ab planes, which makes the value of the orbital 
projection ( (| — 7^ ) (| - rj 7) ) ) in Eq. I|3.12|l very small 
indeed in the physical range of 9 (compare Table [Hj. In 
fact, for the alternating planar (x 2 — z 2 )/(y 2 — z 2 ) orbitals 

one of the operators { (| — t^^ , (| — T j 7 ^) } equals zero, 
and the CT contribution vanishes, so one cannot reduce 
the value of \J a b\ by increasing the AF CT term that 
follows from Ha ■ 

The strong anisotropy of the magnetic exchange in- 
teractions in KCuF 3 is well illustrated in Fig. 01 by the 
ratio \ J a b\/Jc, being close to 0.07 for either the JT OO 




0.00 
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FIG. 4: Ratio of the exchange interactions \J a b\/Jc inKCuF:?, 
given by Eqs. l|IH5}l and for increasing charge transfer 

parameter R, obtained for a few values of n and for the OO 
induced by the JT effect (6 ~ 70°, solid lines). The dashed 
line shows \J a b\/Jc obtained for n — 0.12 and for alternating 
(x 2 - z 2 )/(y 2 - z 2 ) orbitals (6 = 60°). 

(6 ~ 70°), or for the OO suggested by the orbital su- 
perexchange at T < T N (9 = 60°). Note that for R ~ 1 
the ratio \J a b\/Jc does not depend significantly on 9 in 
the interesting range between 9 = 70° and 9 = 60°. 

D. Optical spectral weights for KC11F3 

Now we turn to the optical spectral weights (|2.13|l and 
determine the kinetic energies for the corresponding Hub- 
bard subbands. As discussed in Sec. [H] they originate 
from different multiplet excitations, and depend on the 
OO and on the spin-spin correlations (|2.12|) . Here we 
analyze in detail the spectral weight distribution for po- 
larization along the c axis, where strong exchange in- 
teraction J c controls the spin-spin correlations s c (|2.12() 
which remain finite in a broad temperature regime. 

Knowing that the interchain FM exchange coupling 
J a b is so weak, we describe the temperature variation 
of the spin-spin correlations s c — (Si-Si+i) c employing 
the Jordan- Wigner fermion representation 60 for a ID spin 
chain. One finds for perfect OO at temperature T <T S 

s c = -k(1 + k), (3.27) 

where 

k = — V Icosfcl tanhf-^— V (3.28) 

N ^ 1 1 \2k B TJ v ; 

k 

e k = J c (l + 2/c) I cosfc|. (3.29) 

Here Sk is the ID dispersion of pseudofermions. The ex- 
change interaction J c is constant as long as the orbitals 
remain frozen, and sets the energy scale for the temper- 
ature variation of s c . Eqs. I|3.28|l and (|3.29|) were solved 
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FIG. 5: (Color online) Kinetic energy terms per bond 12.13L 
as obtained for the c axis with the OO 13.141 suggested by 
the JT distortions (9 ~ 70°) in KCuF 3 (in units of J): high- 
spin k[ c ' (solid line) and low-spin K^l (dashed lines), for 
increasing temperature T/J c . Parameters as in Fig. [3] 

self-consistently to obtain s c l|3.27[) as a function of tem- 
perature. In the limit T — ► one finds k = 1/-7T, and 
s c = — (1 + tt)/tt 2 ~ 0.42. This value represents an excel- 
lent analytic approximation to the exact result, 

sr=-(lo g 2-i)~ -0.4431, (3.30) 

obtained for the ID AF Heisenberg chain from the Bethe 
ansatziffi 

The general theory presented in Sec. [n] makes a clear 
prediction concerning the temperature dependence of the 
spectral weights in optical absorption. First of all, a large 
anisotropy between the polarization along the AF c axis 
and the polarization in the (weakly FM) ab plane is ex- 
pected when the AF (FM) spin-spin correlations along 
the c axis (within the ab planes) develop. Indeed, us- 
ing the self-consistent solution of Eqs. (|3.28() and i|3.29[) . 



TABLE III: Kinetic energies of the different Hubbard sub- 
bands (K^), and total kinetic energies (A- (7) ) (in meV), 
as obtained for KCuF3 for two representative orbital states: 
9jt — 70° and #se = 60°. Parameters: J ~ 33.3 meV, 
77 = 0.12. 
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one finds that the kinetic energy K[ (which determines 
the spectral weight of the high-spin excitation at energy 
U — 3Jh) is rather low (Fig. |5J). In contrast, the low-spin 

excitations and contribute with large spectral 
weights in the low temperature regime, reflecting the AF 
correlations along the c axis. 



When the temperature increases, the spin-spin corre- 
lations s c gradually weaken and the kinetic energy redis- 

(c) (c) 

tributes — both low-spin terms K\ ' and K\ ' decrease, 

(c) 

while the high-spin term K\ increases as more high-spin 
excitations are then allowed (Fig. |3J) ■ For the considered 
OO given by 9jt — 70°, the changes of the contributions 
at the two lower energies which correspond to n = 1 
and n = 2 are particularly large between zero and room 
temperature (up to fcsT/J c ~ 0.8), with the increase 

(decrease) of k[ c ^ [K^] by ~ 60 (~ 33) percent of the 
reference value at T = 0. This leads to rather similar- 
values of all three contributions at room temperature. 
This predicted behavior could be verified by future ex- 
periments. 



The temperature variation of the spectral weights for 
ab polarization is more difficult to predict as it involves 
weak s a b spin correlations which develop in the temper- 
ature range T ~ Tjv , and grow with increasing order pa- 
rameter, s a b oc |(S 2 )|, below Tat i^i Assuming the classical 
value for s a b = 0.25 at T — > 0, one finds that the kinetic 
energy would come entirely from the high-spin optical ex- 
citations K[ ab \ while the low-spin excitations would be 
fully suppressed (Table ITTi"|l . Above Tjv the spin system 
is controlled by the dominating AF exchange constant 
J c , and s a b — 0. Even then the high-spin excitations 
at low energy dominate and have large spectral weight 
as a result of the persisting OO. Some decrease of K[ ab ^ 

(c) 

accompanied by the increase of K± ' with increasing tem- 
perature makes the anisotropy between K[ c ^ and K[ ab ^ 
considerably less pronounced, but this anisotropy of the 
spectra in the low energy range remains still close to 3:1 
even at room temperature (Table llllf) . Note also that 
the highest energy excitation for the ab polarization van- 
ishes at 6>se = 60° , and gives a negligible contribution for 
the JT angle 6*jt — 70°, due to the orbital correlations 
within the ab planes. 



Finally, we would like to emphasize again that knowing 
only the exchange constants J a b and J c in KCUF3, one 
is not able to determine all the microscopic parameters 
of the CT model. We emphasize that a better under- 
standing of the properties of KCUF3 can be achieved only 
by combining the results of magnetic and optical experi- 
ments, after the latter experiments have been performed. 
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IV. CUBIC MANGANITE: LaMnO :! 
A. Superexchange model 

Although e g and t 2g electrons behave quite differ- 
ently in LaMnOs and are frequently treated as two 
subsystems) 2 * 3 , the neutron experiments^ which measure 
the spin waves in the A-AF phase below Tn leave no 
doubt that an adequate description of the magnetic prop- 
erties requires a magnetic Hamiltonian of the form given 
by Eq. I|2.11|) . describing superexchange between total 
5 = 2 spins of the Mn 3+ ions. The high-spin 5 E ground 
state at each Mn 3+ ion is stabilized by the large Hund's 
exchange Jh ■ The situation is here more complex than 
either in KCUF3 (Sec. Illlf) or in the t 2g systems discussed 
in the following Sections, however, as the superexchange 
terms between Mn 3+ ions originate from various charge 
excitations d\d 4 =^ d\d^, made either by e g or by t 2g 
electrons, leading to different excited states in the in- 
termediate d 5 configuration on a Mn 2+ (df) ion. Such 
processes determine the U term Hu(d 4 ) defined by Eq. 
(|2.6() , and were analyzed in detail in Ref. 0, and lead to 
the structure of Hu(d 4 ) given below. However, the CT 
processes, H.A(d 4 ), contribute as well and the complete 
model for LaMn03 reads 



n(d 4 ) = nu(d 4 ) + n A (d 4 ). 



(4.1) 



The superexchange constant J is here defined again by 
Eq. (|2.7|l . using an average hopping element along an 
effective (dda) bond, t = t 2 pd jK, where A is an average 
CT excitation energy, introduced below in Eq. (|4.19() . 

First we analyze the structure of the U term for 
LaMnC"3, Hjj{d 4 ), due to excitations involving e g elec- 
trons. The energies of the five possible excited states^ 
(i) the high-spin A\ state(5 = 5/2), and (ii-v) the low- 
spin (5 = 3/2) states: 4 A X , 4 E ( 4 E e , 4 Eg), and 4 A 2 , will 
be parametrized again by the intraorbital Coulomb ele- 
ment U l|2.3l) . and by Hund's exchange Jh between a pair 
of e g electrons, defined in Eq. I|3.2fl . 63 The energies of the 
excited states are given in terms of the Racah parameters 
in Ref. 1391 : in order to parametrize this spectrum by Jh 
we apply an approximate relation AB ~ C which holds 
for the spectroscopic values of the Racah parameters for 
a Mn 2 + (d 5 ) ion^i B = 0.107 eV and C = 0.477 eV. 
Here we use these atomic values as an example of the the- 
ory — using them and Eq. (|3.2[) one finds the excitation 
spectrum: U-3J H , U + SJh/A, £/ + 5J h /4, £/ + 5J h /4, 
and U + 13J ff /4 [Fig. □»]. Unlike J H , the value of U 
is known with less accuracy — hence we shall use it here 
only as a parameter which can be deduced a posteriori 
from the superexchange J which is able to explain the 
experimental values for two exchange constants respon- 
sible for the A-AF phase observed in LaMn03 well below 
the structural transition (here again Tjv €T S ). 

Using the spin algebra (Clebsch-Gordon coefficients), 
and making a rotation of the terms derived for a 
bond (ij) || c to the other two cubic axes a and b, 



one finds five contributions to Tiu(d 4 ) due to differ- 
ent (*l s e g)i(*2g e a)j ^ (*2 ff e g)i(*l fl )j excitations by e g 
electrons li 3 - 



fll 7) =~l(^+6)(±-r 4 ^rj 
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(4.7) 



follow from the above multiplet structure of Mn 2+ (d 5 ) 
ions, and r\ (|3.3|l stands for the Hund's exchange. The 
meaning of the various terms is straightforward: the first 

term describes the high-spin excitations to the 6 Ai 

(7) 

state while the remaining ones, Hn with n = 2, . . . , 5, 
arise due to the low-spin excited states 4 Ai, 4 Eg, 4 E e and 
4 A 2 , respectively. The orbital dependence is given by the 
same operators (13.1111 as in Sec. IIHI Similar to the case 
of the 4 E state for the copper fluoride [see Eqs. I|3.5|) 
and (|3.6p ]. the doubly degenerate 4 E state contributes 
here with two terms characterized by a different orbital 
dependence in Eqs. 14.4(1 and l|4. 5JI . Note that this de- 
generacy would be removed by the cooperative JT effect, 
i.e., the structural phase transition (and associated 00) 
driven by the local JT coupling in combination with the 
elastic lattice forces. The resulting small level splitting 
we neglect here, and so we set r% = r±. 

The superexchange mediated by t 2g electrons re- 
sults from (t| s 4)«(*2g e s)j # (<2 ff e fl)t(^a4)j excitations 
which involve 4 T\ and 4 T 2 configurations at both Mn 
ions: Mn 2+ and Mn 4+ . They give low 5 = 3/2 spins 
of Mn 2+ ions, and this part of the superexchange is AF. 
Using the present units introduced in Eqs. I|2.3|l and 
(13.21) , one finds the excitation energies (not shown in Fig. 
CJ: e( 4 T 1; 4 T 2 ) ~U + 5J ff /4, e( 4 T 2 , 4 T 2 ) ~U + 9J H /4, 
f?( 4 Ti, 4 Ti) ~ U+UJh/A, ande^Ts, 4 ^) ~ U+15J H /4, 
with the first (second) label standing for the configura- 
tion of the Mn 2+ (Mn 4+ ) ion, respectively. In the actual 
derivation each of the excited states, with one t 2g orbital 
being either doubly occupied or empty, has to be pro- 
jected on the respective eigenstates and the spin algebra 
is next used to construct the interacting total 5 = 2 spin 
states. This leads to the final contribution to Tiu which, 
in a good approximation, is orbital independent^ 



H, 



(7) 



Si' St 



(4.8) 
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Here j3 = (t^/t) 2 follows from the difference between the 
effective d—d hopping elements along the a and it bonds, 
and we adopt the Slater-Koster value (3 = 1/9. The 
coefficient r t stands for a superposition of the above ti g 
excitations involved in the ti g superexchange, 



1 



rt = lll 



1 



1 



1 



1 



i 



i 



i 



15„ 
4 



(4.9) 



There is no need to distinguish between the different ex- 
citation energies; all of them arc significantly higher than 
the first low-spin excitation energy for the configuration 
A\, which occurs after an excitation by an e g electron. 

While earlier studies of the superexchange interactions 
in manganites were limited to model Hamiltonians con- 



J 



taining only the U term^i2iiii4 the importance of the 
CT processes was emphasized only recently. 65 For our 
purposes we derived the CT term H.&{d 4 ) by considering 
again excitations by either a or 7r electrons on the bond 
(ij) , leading to two-hole excited states at an intermediate 
oxygen, df{2p {ij) ) 6 dj ^ df (2p {ij) fd]. Unlike in KCuF 3 
with a unique CT excitation, however, in the present case 
a number of different excited states occurs with the exci- 
tation energies depending on the electronic configuration 
of the two intermediate Mn 2+ ions at sites i and j. One 
finds that these various excitations can be parametrized 
by a single parameter R given by Eq. H3.13[) . and the 
excited states on two neighboring transition metal ions 
contribute, as for the U term, both FM and AF terms, 



3; + e) + c 3 (S; -Sj -4 
-cifSi-Sj + e) + c^Si-Sj-A)] ( 

(ij) 



± _ T (7) T (7) 
4 1 3 



Si) 



)} 



where the coefficients c„, with n = 1, • • • , 5, and c t are all determined by 77 and R via rj = rjR: 
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(4.10) 

(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 



The coefficients c n follow from CT excitations by e g elec- 
trons. As in the copper fluoride case (Sec. Illljl . the lowest 
(high-spin) excitation energy will be labeled by A, 



A(%) 



(4.17) 



so the other possible individual (low-spin) excitations at 
each transition metal ion have the energies 

A(%) = A + ^J H , 



A( 4 E) = A + ^-J H , 
A(%) = A + ^J ff , 



(4.18) 



These excitation energies are used here to introduce an 
average CT energy A, 



1 



If- 



10 



(4.19) 



A 26VA( E L4 1 ) A( A A X ) A{ A E) A( 4 A 2 ), 
which serves to define the effective hopping element 



14 



t = tp d /A, and the superexchange energy J 1(2.7(1 in 
a microscopic approach. We emphasize, however, that 
such microscopic parameters as {t p d, A, U, Jh} will not 
be needed here, and only the values of the effective pa- 
rameters {J, 77, R} will decide about the exchange con- 
stants and the optical spectral weights. 

Each coefficient c n (14.11(1 . with n = l,---,5, stands 
for an individual process which contributes with a par- 
ticular orbital dependence due to an intermediate state 
arising in the excitation process, and accompanies either 
a FM or an AF spin factor, depending on whether high- 
spin or low-spin states are involved. As in the U term 
Q4.2 |l -l(4.8 |l . a pair of directional \Q orbitals accompanies 
low-spin excitations, while either high-spin or low-spin 
excited states are allowed when two different orbitals, 
one directional and one orthogonal to it (planar orbital) , 
are occupied at the two Mn 3+ sites. In contrast to the 
U term, also configurations with two planar orbitals |£) 
occupied at sites i and j contribute to Tt A (d ) in Eq. 
(|4.1U|) . These terms are accompanied by the projection 
operator ((± + t^ 7) )(± + rj 7) )) in Eq. jOfy Note that 
in the case of the U term such configurations did not con- 
tribute as the e g electrons were then blocked and could 
not generate any superexchange terms. As the electrons 
from an oxygen 2p a orbital are excited instead to direc- 
tional |C) orbitals at two neighboring Mn 3+ ions, again 
both high-spin (S = 5/2) and low-spin (S = 3/2) excita- 
tions are here possible, giving a still richer structure of 

n A (d 4 ). 

The 00 in LaMn03 is robust and sets in below the 
structural transition at T s ~ 780 KpS& The orbital inter- 
actions present in the superexchange Hamiltonian (|4.1|) 
would induce a G-type 00 £L The observed classical 
ground state, which can again be described by Eq. 1(3.17(1 , 
corresponds instead to C-type 00, as deduced from the 
lattice distortions. Note that in contrast to KCUF3, the 
occupied orbitals refer now to electrons and thus the val- 
ues of the expected orbital angle 9 are > 90° (which 
corresponds to cos 9 < 0) and so are distinctly different 
from the copper fluoride case. The averages of the or- 
bital operators in the orbital ordered state are given in 
Table|nl including the terms ((| + T ( f l) ){\ + rj 7) )) which 
contribute now to the CT part of superexchange. The de- 
pendence on the orbital angle 9 suggests that, similar to 
KCUF3, these new terms are more significant along the c 
axis for the 00 expected in the manganites. 

B. Spin exchange constants and optical intensities 

For a better understanding of the effective exchange 
interactions it is convenient to introduce first the t2 g su- 



perexchange constant Jt which stands for the interaction 
induced by the charge excitations of ti g electrons. When 
the CT terms are included, Jt consists of the two contri- 
butions given in Eqs. ((4.8|) and l(4.1L)|) . 



Jt = Jtu + JtA = gJ0(rt + Ret). (4.20) 



This interaction is frequently called the superexchange 
between the core spins. We emphasize that this term 
is orbital independent and thus isotropic. The coupling 
constant Jt has a similar origin as the e g part of the 

superexchange Je J \ which however is orbital dependent 
and anisotropic. We emphasize that both J t and ji 1 ^ are 
relatively small fractions of J. 

The e g contributions to the effective exchange con- 
stants 1(2.11(1 in LaMn03 depend on the orbital state 
characterized again by Eqs. ((3.17(1 via the averages of 
the orbital operators, 

+ i<* +•»<(! -^)>} 

+ ?(c-c 2 )(i-^>) 

+ ^--K(i+^)(H 7) )>}- <«« 

As the structural transition occurs at relatively high tem- 
perature T s — 780 K, at room temperature (and below 
it) the 00 may be considered to be frozen and specified 
by an angle 9 [see Eqs. I(3.17(l ]. The orbital fluctuations 
are then quenched by the combined effect of the orbital 
part of the e g superexchange in Eqs. I(4.2|) - 1(4.6(1 and the 
orbital interactions induced by the JT effect ji 3 - and the 
spins effectively decouple from the orbitals, leading to the 
effective spin model ((2.11(1 . The magnetic transition then 
takes place within this 00 state, and is driven by the 
magnetic part of the superexchange interactions, which 
follow from Hu(d 4 ) and H A {d 4 ). For the C-type 00, as 
observed in LaMnOs^ 6 . one finds the effective exchange 
constants in Eq. 1(2.11(1 as a superposition of Jt and J^ 7 ' 
after inserting the averages of the orbital operators (see 
Table CU) in Eq. 
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(4.22) 



(4.23) 



Considering /3 = 1/9 to be fixed by the Slater-Koster 
parametrization, a complete set of parameters which de- 
termines J c and J a b comprises: J (|2.7|l . 77 l|3.3|l . R <|3.13|1 . 
and the orbital angle 8 which defines the phase with 00 
by Eqs. (|3.17|) . referring now to the orbitals occupied by 
electrons. 

Equations (|4.22|) and (|4.23|) may be further analyzed 
in two ways: either (i) using an effective model which 
includes only the U superexchange term due to d — d 
transitions, as presented in Ref. [13J and discussed in Ap- 
pendix lAl fi.e.. taking R = which implies TL&(d A ) = 0), 
or (ii) by considering the complete 71(d ) model as given 
by Eq. I|4.1[) . which includes also the CT contributions 
to the superexchange (with R > 0). By a similar anal- 
ysis in Sec. IIIII we have established that the CT terms 
are of essential importance in KCUF3 and should not be 
neglected, as otherwise the strong anisotropy of the ex- 
change constants would remain unexplained. Here the 
situation is qualitatively different — as we show in Ap- 
pendix^ using somewhat modified parameters J and r\ 
one may still reproduce the experimental values of the 
exchange constants, deduced from neutron experiments 
for LaMn03^S within the effective model at R = 0, and 
even interpret successfully the optical spectral 

It is important to realize that the high-spin e g electron 
excitations play a crucial role in stabilizing the observed 
A-AF spin order, as only these processes are able to gen- 
erate FM terms in the superexchange. They compete 
with the remaining AF terms, and the A-AF phase is 
stable only when J a b < and J c > 0. We have ver- 
ified that the terms which contribute to J a b and J c in 
Eqs. H4.22f> and l|4.23[) are all of the same order of mag- 
nitude as all the coefficients {r n ,r t ,c n ,Ct} are of order 
one. Hence, the superexchange energy J <|2.7|l is much 
higher than the actual exchange constants in LaMnO,3, 
i.e., I J a b\ < J and J c < J. 

Next we consider the kinetic energies associated with 
the various optical excitations which determine the opti- 
cal spectral weights by Eq. (|2.13|) . Again, as in the pre- 
viously considered case of KCUF3, the high-spin subband 
at low energy is unique and is accompanied by low-spin 
subbands at higher energies. While the energetic sepa- 
ration between the high-spin and low-spin parts of the 
spectrum ~ 4Jh is large, one may expect that the low- 
spin optical excitations might be difficult to distinguish 
experimentally from each other. As we will see below by 
analyzing the actual parameters of LaMnOs, the low-spin 
excitations overlap with the d — p CT excitations, and so 
such a separation is indeed impossible — thus we analyze 



here only the global high-energy spectral weight due to 
the optical excitations on the transition metal ions, ex- 
pressed by the total kinetic energy Kf^ for all (e g and 
t2 g ) low-spin terms, and compare it with the high-spin 

one, given by k[ 7 \ Using the manganite model 14. 1|) 
one finds for polarization in the ab plane 
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and for polarization along the c axis 
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The optical spectral weights given by Eqs. (|4.24|) - (|4.27l) 

are determined by J, 77, and the orbital angle 9. Note that 
the leading term in the low-spin part comes from the e g 
optical excitations, while the t2 g excitations contribute 
only with a relatively small weight oc (3. 



C. Magnetic interactions in LaMnC>3 

It is a challenge for the present theoretical model to 
describe both the magnetic exchange constants 62 and the 
anisotropic optical spectral weights 28 using only a few 
effective parameters {J, 77, R} and the orbital angle 9. We 
shall proceed now somewhat differently than in Sec. IIIII 
and analyze the experimental data using primarily these 
effective parameters, while we will discuss afterwards how 
they relate to the expectations based on the values of the 
microscopic parameters found in the literature. 

The experimental values of the exchange constants^ 
J a b = —1.66 meV and J c = 1.16 meV, impose rather 
severe constraints on the microscopic parameters and on 
the possible OO in LaMnOs. The AF interaction J c is 
quite sensitive to the type of occupied orbitals l|3.17|l . 
and increases with increasing amplitude of \z) orbitals 
in the ground state, i.e., with decreasing orbital angle 
9. Simultaneously, the FM interaction — J a b is enhanced. 
Already with the effective model (at R = 0) it is not 
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straightforward to determine the parameters J and rj, as 
we discuss in Appendix [X] This model is in fact quite 
successful, and a reasonable agreement with experiment 
could be obtained both for the magnetic exchange con- 
stants and for the optical spectral weights, taking the 
experimental excitation spectrumtSS Here we will inves- 
tigate to what extent this effective model gives robust 
results and whether the CT processes could play an im- 
portant role in LaMnC>3. 

By analyzing the CT terms in the superexchange [com- 
pare Eqs. I)4.22[l . I|4.23jl ] one concludes that these con- 
tributions are predominantly AF. Therefore, it might be 
expected that a higher value of 77 than 0.16 used in Ap- 
pendix ^ would rather be consistent with the magnetic 
experiments. Increasing 77 gives an increased coefficient 
7*1, so not only the FM term in J 7 is then enhanced, but 
also the optical spectral weight K f b which corresponds to 
the high-spin transition. Simultaneously, the angle 8 is 
somewhat increased, but the dependence of the spectral 
weight K i b on the angle 8 is so weak that the higher value 
of 77 dominates and a somewhat lower value of J than 170 
meV used in Appendix [X] has to be chosen. Altogether, 
these considerations led us to selecting J = 150 meV and 
r\ = 0.181 as representative parameters for which we show 
below that a consistent explanation of both magnetic and 
optical data is possible. 

After these two parameter values have been fixed, it 
is of interest to investigate the dependence of the effec- 
tive exchange interactions on the CT parameter R and 
on the orbital angle 8. As in the KCUF3 case, one finds 
a stronger increase of J c with increasing R, while these 
terms are weaker and lead to nonmonotonic changes for 
Jab (Fig. EI). First of all, with the present value of 
?; = 0.181, at R = the AF interaction J c is close to zero 
for angles 8 ~ 100°, while the FM interaction J ao < —2 
meV is somewhat too strong. With increasing R one 
finds that J c increases, while the FM interaction J a ;, ini- 
tially becomes weaker when R increases up to R = 0.3 
and the term oc c± dominates the CT contribution to J a b 
[see Eq. (|4.22|l ]. At higher values of R, however, the FM 
contributions due to C3 — C2 < and C5 — C4 < start 
to dominate, and J a b decreases with increasing R, par- 
ticularly for small values of 8 < 90°. One finds that the 
experimental values of both exchange constants are well 
reproduced for R — 0.6 and at the orbital angle 8 = 102°. 

Although this fit is not unique, one has to realize that 
the experimental constraints imposed on the parameters 
are indeed rather severe — as we show below the present 
parameters give a very reasonable and consistent inter- 
pretation of the experimental results for LaMnC>3. For 
the above parameters the FM and AF terms to J c almost 
compensate each other in the U term, and a consider- 
able AF interaction along the c axis arises mainly due to 
the CT contributions (Fig. 0). This qualitatively agrees 
with the situation found in KCUF3, where the CT term 
was of crucial importance and helped to explain the ob- 
served large anisotropy between the values of exchange 
constants. Also the CT term which contributes to J a b 
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FIG. 6: (Color online) Superexchange constants for LaMnOs: 
(a) J c along the c axis, and (b) J a b in the ab plane, as func- 
tions of orbital angle 9 which defines the C-type OO, see Eqs. 
13.171 . Dashed lines for increasing R = 0, 0.15 and 0.30 from 
bottom to top; solid lines for R — 0.6. Note that the depen- 
dence of J a b on R is nonmonotonic. Horizontal dotted lines 
indicate the experimental values of Ref. |62| for LaMnOa; di- 
amond and circle show the experimental values of J c and J a b 
for 6 = 102°. Parameters: J = 150 meV, r/ = 0.181. 



is AF and increases with increasing angle 8, while the U 
term is FM but weakens at the same time. This results 
in a quite fast dependence of the FM interaction J a b on 
8 (Fig. 0. 

One thus recognizes a similar dependence of the ex- 
change constants — J a b and J c on the orbital angle 8 
(Fig. Q) as that found before with the U term alone 
(see Fig. ^]in Appendix The CT terms have mainly 
two consequences: (i) a large AF contribution to the in- 
teraction along the c axis J c , and (ii) an increase of the 
orbital angle well above 8 = 90°. These trends are robust 
in the realistic parameter range. Therefore, one expects 
that the occupied e g orbitals approach the frequently as- 
sumed alternating directional orbitals in the ab planes 
(3x 2 - r 2 )/(3y 2 - r 2 ) with 8 = 120°, but cannot quite 
reach them. Indeed, we have verified that orbitals with 
8 ~ 120° are excluded by the present calculations, as 
then the FM interaction within the ab planes changes 
sign and becomes weakly AF. Thus, the mechanism for 
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pound considered in this paper, due to the fact that the 
d shell is half-filled in the Mn 2+ excited states^ This 
leads to relatively large FM contributions, even when the 
orbitals are not so close to the ideal alternation of direc- 
tional and planar states (as found along the c axis) , which 
would maximize the averages of the orbital operators, 
(j — T i T j)i that control the weight of the high-spin 
terms, see Table ITT1 This result is remarkable but again 
qualitatively the same as found in the effective model 
of Appendix IX! However, quantitatively the ti g term is 
here somewhat stronger, as Jt is now increased by ~ 35% 
over its value J t = 1.7 meV deduced from the effective 
d — d model with U terms only. As a common feature 
one finds that J t ~ 2 meV, so we emphasize that the 
superexchange promoted by ti g electrons is quite weak 
and is characterized by a small value of Jt ~ 4 x 10~ 3 i 
(or ~ O.Olt for the superexchange between S = 3/2 core 
spin a 2 ' 71 ) which provides an important constraint for the 
realistic models of manganites. 



D. Optical spectral weights in LaMnOa 



FIG. 7: (Color online) Superexchange interactions — J a b 
14. 2211 and J c 14.231 1 as functions of the orbital angle 6, ob- 
tained within the superexchange model l|4.1|l for LaMnOs 
(solid lines) . Contributions due to the U term and due to the 
CT term to both exchange constants {—J a b and J c ) are shown 
by dashed and dashed-dotted lines, respectively. Isotropic and 
orbital independent tig terms superexchange terms Jtu and 
JtA 14.201 1 are: 1.45 and 0.83 meV. Experimental values 62 
of —Jab and J c (indicated by diamond and circle) are well 
reproduced at angle 8 — 102°. Parameters J = 150 meV, 
■q = 0.181, R = 0.6. 



the observed ^4-AF phase is lost, and one has to conclude 
that angles 9 > 9jt are excluded. 

Indeed, we have found that the orbital angle 9 ~ 102° 
reproduces well the experimental data for both exchange 
constants (Fig. 0), and is thus somewhat smaller than 
the angle 0jt = 108° deduced from the lattice distor- 
tions in LaMnOs^ This can be seen as a compromise 
between the orbital interactions involving the lattice and 
the purely electronic superexchange orbital interactions, 
so it is reasonable to expect that 6 < 9jt- 

Finally, we emphasize that the e g electron excitations, 
contributing both to the U and to the CT processes in 
the superexchange, are FM for all cubic directions^ 
and only due to a substantial tig term which follows 
from low-spin (AF) excitations, J t = 2.28 meV, the ex- 
change interaction along the c axis changes its sign and 
becomes AF. Alltogether, the present analysis shows that 
the t^g superexchange plays a decisive role in stabilizing 
the observed A-AF spin order — without it already the 
undoped LaMnOa would be a ferromagnet. This pecu- 
liar situation follows from the large splitting between the 
high-spin and low-spin excitations ~ 3 eV in LaMnO,3, 
which is larger than in any other transition metal corn- 



Consider now the temperature dependence of the op- 
tical spectral weights (12.15(1 . As the orbital dynamics is 
quenched up to room temperature T ~ 300 K, it suf- 
fices to consider the temperature dependence of the in- 
tersite spin-spin correlations l|2.12|l . We derived these by 
employing the so-called Oguchi method 7 ^ (see Appendix 
lB)l . which is expected to give rather realistic values of 
spin correlations functions for the large spins S = 2 in 
LaMnOs in the entire temperature range. Thereby we 
solved exactly the spin-spin correlations {s a b,s c } on a 
single bond (ij) coupled to neighboring spins by mean- 
field (MF) terms, proportional to the order parameter 
(S z ). A realistic estimate of the magnetic transition can 
be obtained by reducing the MF result obtained for S = 2 
spins by an empirical factor ~ 0. 71 M> Using the exchange 
interactions obtained with the present parameters, one 
finds Tn ~ 143 K which reasonably agrees with the ex- 
perimental value of T™ p = 136 Km£ The calculations of 
spin-spin correlations are straightforward, and we sum- 
marize them in Appendix iBl Both correlation functions 
s a b and s c change fast close to Tjv, reflecting the temper- 
ature dependence of the Brillouin function which deter- 
imes (S z ), and remain finite at T> T/y. 

The large splitting between the high-spin and low-spin 
excitations makes it possible to separate the high-spin ex- 
citations from the remaining ones in the optical spectra, 
and to observe the temperature dependence of its inten- 
sity for both polarizations. As in the effective modelj22iZi 
the present theory predicts that the low-energy optical 
intensity exhibits a rather strong anisotropy between the 
ab and c directions. It is particularly pronounced and 
close to 10:1 at low temperatures when the spin correla- 
tions \s a b\ and s c are maximal (Fig. |SJ). Unfortunately, 
this anisotropy at T — > in only weakly dependent on 
the orbital angle 9, so it cannot help to establish the type 
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FIG. 8: (Color online) Anisotropy of the low-energy spec- 
tral weights in LaMnOa, given by the ratio K[ ab) / K[ c \ for 
increasing temperature T, as obtained for: 8 = 102° (solid 
line), 6 = 108° (dashed line), = 98° (dashed-dotted line). 
The arrows indicate the Neel temperature Tn derived from 
the exchange constants in each case (see text). Parameters: 
J = 150 meV, rj = 0.181. 



of OO realized in the ground state of LaMnOa, and its 
possible changes with increasing temperature. In fact, 
when the parameters are fixed and only the orbital angle 
9 is changed, a different value of the Neel temperature 
follows from the modified exchange constants, being the 
main reason behind the somewhat different temperature 

dependence of the ratio K^/K^ (Fig. |SJ). 

One finds a very satisfactory agreement between the 
present theory and the experimental results of Ref. l28l 
as shown in Fig. [5] We emphasize, however, that 
the temperature variation of the optical spectral weights 
could also be reproduced within the effective model at 
R = O.Siii showing that the CT terms lead only to minor 
quantitative modifications. Note also that at this stage 
no fit is made, i.e., the kinetic energies (|2.13l) which stand 
for the optical sum rule are calculated using the same pa- 
rameters as found above to reproduce the exchange con- 
stants in Fig. Therefore, such a good agreement with 
experiment demonstrates that indeed the superexchange 
interactions describe the spectral intensities in the optical 
transitions. We note, however, that the anisotropy in the 
range of T > TV is somewhat larger in experiment which 
might be due to either some inaccuracy of the Oguchi 
method or due to the experimental resolution. 

The distribution of the optical intensities and 
their changes between the low (T ~ 0) and the high tem- 
perature (T = 300 K) regime are summarized in Table 
IIVI At temperature T = one finds that for the ab polar- 
ization the entire spectral intensity originates from high- 
spin excitations. This result follows from the classical 
value of the spin-spin correlation function s a f, = — 1 pre- 
dicted by the Oguchi method. As quantum fluctuations 
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FIG. 9: (Color online) Kinetic energies per bond (in meV) 
as functions of temperature T, obtained for LaMnOs for the 
high-spin (solid lines, n = 1) and for the low-spin (dashed 
lines, n > 1) excitations at the Mn ions for polarization: (a) 



in the ab planes, K$? b \ and (b) along the c axis, K^ c> . The 
total kinetic energies are given by the long dashed-dotted 
lines, while the vertical dotted lines indicate the value of Tn 
derived within the present model (see text). Filled circles 
show the experimental intensities— at low energy (n = 1). 
Parameters: J = 150 meV, i] ~ 0.181, R = 0.6, and ~ 102°. 



(c) 



in the ^4-AF phase are small^ the present result is nearly 
accurate. When the temperature increases, one finds con- 
siderable transfer of optical spectral wei ghts between low 
and high energies, discussed also in Ref. |2£|, with almost 
constant total intensities K^ 1 ' for both 7 = ab and c (see 
also Fig. The optical weights obtained for the JT an- 
gle #jt = 108° are similar to those for 6 — 102°, showing 
again that the optical spectroscopy is almost insensitive 
to small changes of the orbital angle. In contrast, for 6>jt 
the exchange constants (and hence also the estimated 
value of Tat) are too low. 

While the effective parameters {J, rj,R} used in this 
Section give a very satisfactory description of both mag- 
netic and optical properties of LaMn03, the values of 
the microscopic parameters, such as the Coulomb inter- 
action on the Mn ions U, and on the oxygen ions U p , 
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TABLE IV: Exchange constants J 7 [see Eqs. 14.221 and 
H, kinetic energies for high-spin (K^) and for low-spin 
(iff,g ) excitations [see Eqs. I4.24l - 14.27l l. and total kinetic 
energies K[2l (all in meV), for a bond within an ab plane 
(7 = ab) and along a c axis (7 = c), as obtained for the OO 
given by the orbital angle 6 ~ 102° which reproduces the ex- 
perimental exchange constants, 63 : and for the angle suggested 
by the lattice distortions, 9jt- Parameters: J = 150 meV, 
7/ ~ 0.181, R = 0.6. 
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the Hund's exchange J#, the charge transfer parameter 
A, and the d — p hopping element t p d, are not uniquely 
determined. One could attempt to fix these parameters 
using the atomic value of Hund's exchange Jh — 0.9 eV. 
With T] = 0.181 it leads to U ~ 5.0 eV, quite close to 
other estimates, 43 while the value of R = 0.6 suggests 
then that, taking the commonly accepted 4 ^*^ value of 
the CT energy A = 5.0 eV, the oxygen Coulomb ele- 
ment is large, U p ~ 7 eV. This value of U p is larger 
than usually obtained U p *~ 4 eV for oxygen ions, as for 
instance by analyzing the parameters of the three-band 
model for Cu02 planes. 5 — We note, however, that the op- 
tical data^ suggest a somewhat reduced value of Hund's 
exchange Jh ~ 0.7 eV (using the present units), so it 
could be that the local exchange interactions are some- 
what screened in reality. 

Fortunately, when the optical data are available, also 
the position of the low-energy excitation is known, and 
this may serve as an additional experimental constraint 
for the parameters^! This excitation is found at about 
2.0 eV^ indicating that U - 3J H ^ 2.0 eV. With this 
constraint one finds, using again 77 = 0.181, that U ~ 3.7 
eV and Jh — 0.67 eV. These parameters give the low- 
spin A A\ and 4 E excitations close to 4.5 eV, in agreement 
with experiment. 28 So the above values of the microscopic 
parameters {U, Jh} appear to be consistent both with 
the present value of 77 and with the spectra. Further- 
more, for these empirical parameters one finds R = 0.6 
either with A = 5.0 and U p ~ 2.3 eV, or with A = 4.2 
and Up ~ 4.0 eV. These values, particularly the second 
set, are perfectly acceptable and in the usually consid- 
ered rangei 54 ' 64 Taking the above value of U ~ 3.7 and 
J = 150 meV, one finds t = 0.37 eV, a somewhat lower 
value than that which follows from the effective model of 
Appendix [S] Altogether, these results indicate, contrary 



to what is frequently assumed* 4 ^* that the local ex- 
change interactions are somewhat screened in reality by 
covalency effects, and that at the same time the screen- 
ing of the intraorbital Coulomb interaction U is stronger 
than estimated beforei 4 ^^ 



V. CUBIC TITANITES 
A. Spin-orbital superexchange model 

Perovskite titanates, LaTiOs and YTi03, are intrigu- 
ing examples of Mott insulators with orbital degrees of 
freedom due to ti g electrons: in the ground state the 
Ti 3+ ions are in the t\ g configuration. In an ideal per- 
ovskite structure the ti g orbitals are degenerate, but lat- 
tice distortions may contribute to the magnetic ground 
stato 76 i 77 i 78 i 79 and to the Mott transition^ — here we 
do not intend to discuss this controversial issue. 

In an ideal cubic system each ti g orbital is perpendic- 
ular to a single cubic axis, for instance the \xy) orbital 
lies in the ab plane and is perpendicular to the c axis. It 
is therefore convenient to introduce the following short 
hand notation for the orbital degree of freedom^ 

\a) = \yz), \b) = \zx), \c) = \xy). (5.1) 

The labels 7 = a, 6, c thus refer to the cubic axes perpen- 
dicular to the planes of the respective orbitals. 

The superexchange spin-orbital model (|2.8|) in cubic ti- 
tanates couples 5=1/2 spins and the orbital ti g degrees 
of freedom at nearest neighbor Ti 3+ ions. Due to large U 
the electron densities satisfy thereby the local constraint 
at each site i, 

n la + n ib + n ic = 1. (5.2) 

In titanates there is no need to consider CT processes, as 
these systems are Mott-Hubbard insulators 1 " and no qual- 
itatively new effects apart from some negligible renormal- 
ization of the effective parameters {J, 77} could arise from 
CT excitations. This simplifies our considerations, so we 
analyze the superexchange in the leading order of per- 
turbation theory, given by the contributions which result 
from virtual excitations between the neighboring Ti 3+ 
ions, {t\ g )i{t\ g )j =^ (^g)i(i§ff)j- These charge excitations 
involve the Coulomb interactions in the d 2 configuration 
of a Ti 2+ ion, parametrized as before by the intraorbital 
Coulomb element U, and by Hund's exchange element Jh 
for a pair of ti g electrons, defined as follows (see Table 

Jh = 3B + C. (5.3) 

The charge excitations lead to one of four different ex- 
cited states^ 9 - shown in Fig. db): the high-spin 3 Ti state 
at energy U — 3Jh, and three low-spin states — degen- 
erate 1 T2 and 1 E states at energy U — Jh, and an 1 A\ 
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state at energy U + 2Jh- As before, the excitation ener- 
gies are parameterized by ij, defined as in Eq. (|3.3|) . and 
we introduce the coefficients 
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One finds the following compact expressions for the su- 
perexchange Huid 1 ) as given by Eq. (|2.6|) : 15,81 
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As in Sees. Ill and IV the orbital (pseudospin) op- 
erators {A^J\b^7 , nff} depend on the direction of 
the (ij) bond. Their form follows from the simple 
observationiS^ that only two i 2g orbitals (flavors) are 
active along each cubic axis, e.g. if 7 = c, the active 
orbitals are a and b [see Eq. QS.lfl ]. and they give two 
components of the pseudospin T = 1/2 operator fi. Here 

the operators {A™ , B^l'} describe the interactions be- 
tween the active orbitals along the particular bond, which 
include the quantum effects due to their fluctuations, and 



take either the form of a scalar product f 
a similar expression, 



fj , or lead to 



T- T ■ 

'» 3 



T y T V 
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(5.11) 



which involves double excitations due to t, + t, + and t- t 
terms (as in the U < Hubbard model). The interac- 
tions along axis 7 are tuned by the number of electrons 



,(7) 



1 



which is fixed 



occupying active orbitals, 
by the number of electrons in the inactive orbital rii 7 , 
because of the constraint ^5.2(1 . 



B. Spin exchange constants and optical intensities 

The exchange constant for a bond (ij) along axis 7 
is obtained from Eqs. H5.5[) - l|5.7[) by averaging over the 
orbital states at both sites i and j, 



\j{n+r,){^)- l -J(r 2 -n){B^) 
jJ(ri-r 2 )(r4 7) +nJ 7) ). (5.12) 



The cubic titanates are known to have particularly 
pronounced quantum spin-orbital fluctuations, and 
their proper treatment requires a rather sophisticated 
approachii^iSi Here we shall ignore this complex quan- 
tum problem, and shall illustrate the general theory by 
extracting the magnetic exchange constants from Eq. 
()5.12(l . and the optical spectral weights Ij2.15|l . using an 
ansatz for the 00 in the ground state, in analogy to the 
approach employed in Sees. IIIII and lTvl ior the more clas- 
sical e g systems. 

In general a classical orbital state in the titanates with 
GdFe03-type lattice structure can be parametrized as 
follows, 

IV>i> = a\a)+(3\b)+j\c), 
|V> 2 ) = (3\a)+a\b) + 7 \c), 
^3) = a\a)+f3\b)-l\c), 
^4) = 0\a)+a\b)—y\c), (5.13) 



with real coefficients and normalized (a 2 + (3 2 + j 2 



1) 



wavefunctions at each site. The occupied orbitals \ipi) 
refer to four sublattices (i = 1, •••,4), with Tii(OOO), 
Ti 2 (100), Ti 3 (001), Ti 4 (101) positionsi&S The minus 
signs in \tps) and |V>4) reflect a mirror symmetry present 
in the GdFe03 structure. Note that this state resembles 
G-type 00, and is thus different from the C-type 00 
encountered for e g orbitals [due to the change of sign for 
the \c) orbitals in Eqs. (|5.13l) along the c axis]. 

Using the ansatz l|5.13|l one finds after a straightfor- 
ward calculation the exchange constants, 



Jc 



Jab — 



1 



-J 



Oi + r 2 r 3 )(l- 7 2 ) 2 
"(ri-r 2 )(l- 7 2 ) , 

2(ri+r 2 r 3 )(7 2 + a£) 2 
-{r% -r 2 )(l+ 7 2 ) . 



(5.14) 



(5.15) 



They are determined once again by: (i) the superex- 
change parameter J given by Eq. I|2.7|l . with t standing 
now for the effective (ddir) hopping element, (ii) Hund's 
exchange element rj (13.31) . 84 ° and (iii) the orbital state 
(|5.13|) . specified in the present case by the coefficients 
{a,/3,7}. 

Following the general theory of Sec. |Hj the optical ex- 
citations corresponding to the high-spin Hubbard band 
at energy U — 3Jh (n = 1), and to the two low-spin Hub- 
bard bands at U — Jh (n — 2) and U + 2Jh (n = 3) of 
Fig. njb), have total intensities given by the respective 

kinetic energies see Eq. (|2.15|l . Using the classical 

wavefunctions given by Eqs. (|5.13|) . one finds the follow- 
ing optical spectral weights: 
— for polarization along the c axis, 



K 



(<=) 



K, 



^ l7 2 (l-7 2 )(| + sc), 
^ 2 (l-7 2 )(l^7 2 )( 



(5.16) 
(5.17) 
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K ( 3 C) = ^Jr 3 (l- 7 2 ) 2 (i- Sc ), (5.18) 
— for polarization in the ab plane, 

= Jri [^(l+7 2 ) - (7 2 + «/3) 2 ] (J , (5.19) 

= ^ 2 [i(l+7 2 ) + \{l 2 + ^) 2 ] (j - Sab ),(5.20) 
^ ab) = ^r 3 ( 7 2 +a/3) 2 (i- Safc ). (5.21) 

The kinetic energies ifi 7 ^ depend on the same parameters 
{J, 77} as the exchange constants (|5.14|) and (|5.15l) . on the 
orbital state l|5.13[) via the coefficients {a,/3, 7}, and on 
the spin-spin correlations l|2.12ll . 



C. Application to LaTi0 3 and YTi0 3 

It is now straightforward to investigate the dependence 
of the magnetic and the optical properties in the four- 
sublattice classical state (|5.13[1 on the effective parame- 
ters {J, 77} and on the type of 00 given by the coefficients 
{a, (3, 7}. A priori, the average electron density in the 
|c) orbitals, n c — {rii C ), is different from the densities in 
the other two orbitals (n a and rib), and the cubic sym- 
metry of the expectation values J 7 and K% is explicitly 
broken by the 00 given by Eq. 15.13fl . unless all the 
orbital amplitudes are equal, 



0.5 



a = /3 = 7 = 7f 



(5.22) 



as argued in Refs. \7&7K First we consider this isotropic 
state and evaluate the exchange constants J 7 for increas- 
ing Hund's parameter 77. One finds then AF superex- 
change J 7 ~ 0.4 J which decreases with increasing 77 [Fig. 
llDf a)]. but the interactions remain AF in the physically 
relevant range of 77 < 0.28. Note that this classical con- 
sideration seriously overestimates the actual exchange in- 
teraction as one expects instead J 7 ~ 0.16J when quan- 
tum effects are included^ Nevertheless, having no in- 
formation about the optical excitations, we give here an 
example of estimating the exchange constants and opti- 
cal spectral weights using the spectroscopic values^ for 
Hund's exchange Jh — 0.59, and for Coulomb intraor- 
bital interaction of U = 4.8 eV, which give 77 = 0.123. 
Assuming now a hopping parameter t — 0.2 eV, one finds 
J ~ 33 meV which gives J 7 ~ 14 meV. While one might 
expect that the effective U is somewhat reduced as in the 
case of LaMnOa, and the value of 77 would then be larger, 
the present crude estimate is accidentally quite close to 
the experimental value of J cxp = 15.5 meVS 

In the AF state, realized in LaTiOa, all three Hubbard 
subbands contribute in the optical spectroscopy, and tak- 
ing the Neel state with s a i, = s c = —0.25 one finds the 
highest spectral weight at T = for n = 2 [Fig. fHfb)]. 
For the realistic value 77 ~ 0.123 the spectral weight of 
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FIG. 10: (Color online) Magnetic and optical properties of 
LaTiOs, as obtained for the classical wavefunctions 15. 1311 
with isotropic orbital amplitudes 15.221 , for increasing Hund's 
exchange rj = Jh/U: (a) exchange interactions J c = J a j,, and 
(b) kinetic energy terms K„'^ at T = 0: high-spin [n = 1, 
solid line), and low-spin intermediate-energy (n = 2, dashed 
line) and high-energy excitation (n = 3, long-dashed line). 
The total kinetic energy is shown by the dashed-dotted 
line. 



the lowest-energy (high-spin) Hubbard subband is 
relatively weak, and is similar to that of the n = 3 low- 
spin excitation at the highest energy [Fig. HOT b)]. Note, 
however, that in early optical experiments the intensity 
was found to be practically independent of energy so 
different excitations might be difficult to separate from 
each other. 

Let us verify now whether the wavefunctions as given 
by Eqs. (|5.13l) could also lead to isotropic AF states for 
other choices of orbital amplitudes than in Eq. Q5.22fl . 
Thus, we considered equal amplitudes a and 7 in the 
states parametrized by 5, 

a = 7 = (l-*)-L + a-l (5.23) 



with < 
/?=(!- 



S < 
2 



1. The normalization condition gives 
7 2 ) 1 / 2 which vanishes at S = 1. As 
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FIG. 11: Exchange interactions J c and J a b, Eqs. 15.121 (solid 
and dashed line) for the titanite model for different types of 
orbital order given by Eqs. 15.1311 . The coefficients in the 
wavefunctions l|5.130 vary between the isotropic state a = 
(9 = 7= 1/V3 and the state with linear combinations of 
\a)/\b) and |c) orbitals, with the coefficients given by: a = 7 = 
l/v / 3 + 5(l/v / 2-l/\/3) and ,3 = (l-a 2 -7 2 ) 1/2 . Parameters: 
■q = 0.123. 

expected, the exchange interactions are anisotropic for 
small S > (Fig. The increasing occupancy of the 

c orbitals with increasing 8 results in a somewhat en- 
hanced exchange interaction J ao near the isotropic state 
(for 5 < 0.5), while the interaction J c decreases almost 
linearly^ At larger orbital anisotropy the interaction J ao 
decreases as well. It is remarkable that the AF interac- 
tions become again isotropic close to a = 7 = 1/V2 
which gives a large charge density n c ~ 0.5 in the |c) 
orbitals, but this feature is not generic and will likely be 
modified by quantum effects. Note also that the state 
at (3 = in not a FM state, contrary to some earlier 
suggestions^ 

Finally we remark that the temperature dependence 
of the optical spectral weights can again be studied in 
a similar way as for LaMn03, by solving the spin-spin 
correlation functions at increasing temperature using the 
method of Appendix iBl for S = 1/2 spins. Having no 
experimental data neither for LaTiC>3 nor for YTiC>3, 
we shall limit ourselves to making predictions concerning 
the overall change of the spectral weights between low 
temperature T ~ where the spin-spin correlations are 
maximal, and room temperature T 3> Tjv (or T 3> Tc), 
where the spin-spin correlations may be neglected. We 
take rigid classical wave functions with OO (|5.13(l and 
discuss first AF states. Apart from the isotropic AF 
phase H5.22fl . called state A, we consider in Table IVl also 
somewhat modified orbital amplitudes, as proposed by 
Itoh et alM 

a = 0.690, (9= 0.452, 7 = 0.565, (5.24) 
called here state B. One finds that the exchange inter- 



actions J c and J ao are then rather anisotropic (see Table 
|VJ, so we believe that this state can be experimentally 
excluded on the basis of the neutron data^ 

Due to the AF order in the ground state for the wave- 
functions A and B, the lowest energy high-spin excita- 
tions are suppressed at low temperature, but we predict 
that their weight should increase by about 50% above the 
magnetic transition (Table IV|l . An even stronger tem- 
perature dependence is predicted for the low-spin part 
of the spectra with its intensity decreasing by a factor 
2 when the magnetic order is lost. Thus, the spectral 
weight is shifted from the high-energy to the low-energy 
part of the optical spectrum with increasing temperature, 
and the total spectral weight decreases, both features be- 
ing typical of AF bonds. An experimental confirmation 
of this trend would also allow one to identify the en- 
ergy splitting between the high-spin and low-spin parts 
of the upper Hubbard band, and to determine the value 
of Hund's exchange Jh from experiment. Although the 
anisotropy in spectral weights found in the state (|5.24() 
proposed by Itoh et alJ— might be too weak to be seen 
experimentally, the predicted anisotropy in the exchange 
interactions is certainly detectable. 

It is quite remarkable that an isotropic FM ground 
state of YTi03 can also be described classically by the 
OO state with four sublattices, but the phase factors in 
the wavefunctions have to be then selected differently: 81 

l&> = ^(|o> + l&> + |c», 
\<h) = -*=(-|a)-|&> + |c», 
\<h) = ^(-|o) + |6)-|c)), 

TABLE V: Exchange constants J 7 , kinetic energies for 
Hubbard subbands K„ and total kinetic energies (all 
in units of J) for two orbital ordered states suggested for 
LaTiC>3: state A 15.221 . and state B 15.241 1 reported in Ref. 
IH, and for the OO state suggested for YTi0 3 in Ref. 
Parameter: r\ = 0.123. 







LaTi0 3 




YTiOa 




state A <5"23l 


state B jS^H 


state F JS^lpt 




T = 


T>Tjv 


T = 




T = 


T>T C 


Jc 


0.407 


0.428 


-0.148 


Jab 


0.407 


0.351 


-0.148 




0.176 


0.264 


0.172 


0.258 


1.057 


0.792 




0.465 


0.232 


0.476 


0.238 


0.000 


0.190 




0.118 


0.059 


0.124 


0.062 


0.000 


0.000 


K (c) 


0.759 


0.555 


0.772 


0.558 


1.057 


0.982 


K (ab, 


0.176 


0.264 


0.207 


0.311 


1.057 


0.792 


K (ab) 


0.465 


0.232 


0.451 


0.226 


0.000 


0.190 


K (ab) 


0.118 


0.059 


0.107 


0.053 


0.000 


0.000 




0.759 


0.555 


0.765 


0.590 


1.057 


0.982 



23 



|0 4 > = ^=(|a) - \b) - |c». (5.25) 

We call this state the orbital F state. Neglecting 
again orbital (quantum) fluctuations, one finds (A^ ) = 

(B^) = in this state, while raj 75 = |. At finite 77 
this leads indeed to an isotropic FM state with exchange 
constants, 

7™ = ^ = --./^, (5-26) 

in lowest order oc 77^ so this type of 00 leads to a 
markedly different exchange constant from that given by 
Eq. (|5.14() . Assuming again the same parameters as 
above one finds J^ M ~ —5 meV, which is somewhat 
higher than the experimental value J° xp ~ — 3 meV. 90 
However, in view of the simplicity of the present consid- 
erations this agreement can be regarded as satisfactory. 

Consider now the spectral weights in the optical spec- 
troscopy. One finds rather simple expressions for the 
(isotropic) optical spectral weights for the F state: 

K ™ = 5 Jri (l + s 0' (5 ' 27) 

= 0. (5.29) 

Thus, the optical spectrum of YTi03 would be quite dif- 
ferent from that of LaTi0 3 (see Table EJ- For the FM 
state at T = (with s 7 = j) only the high-spin spec- 
tral weight = I Jr\ contributes, and for the present 

(7) 

parameters K\ ~ J. Thus, future optical experiments 
on YTi03 combined with magnetic experiments could 
help to determine the superexchange constant J. Fur- 
thermore, for the FM F orbital state (|5.25|l one finds at 
high temperature a moderate reduction of the low-energy 
spectral weight by ~ 25%, and a corresponding increase 
for the higher energy excitation n = 2 [the third excita- 
tion does not contribute as long as the 00 remains the 
same as in Eqs. Q5.25[l ]. The energy difference between 
these excitations is 2Jh (Fig. QJ, so they might serve to 
determine this parameter from the optical spectroscopy. 

The above results show once again that the magnetic 
interactions and the optical spectral weights depend in a 
crucial way on the underlying orbital state. Once the 00 
in F state (|5.25|) has been fixed, the exchange constants 
J™ (|5.26[) were completely determined by the superex- 
change J H2.7fl . and by Hund's exchange parameter 77. 
For the optical spectral weig one 
needs in addition the spin-spin correlation function s 1 . 
Most importantly, a strong temperature dependence of 
the optical spectra follows here from a classical orbital 
picture, and this should be a crucial experimental test of 
the validity of the present approach. 



VI. CUBIC VANADATES 
A. Spin-orbital superexchange model 

The last example of spin-orbital physics we want to 
give here are the cubic vanadates, in which the orbital 
degrees of freedom originate from the t 2g configuration 
of the V 3+ ions present in the ground state. Here our 
goal will be to answer the question to what extent these 
systems could be described by decoupling spin and or- 
bital degrees of freedom and assuming classical states 
with 00. The 00 in the cubic vanadates is complemen- 
tary to the spin order — the C-type AF phase is found 
for G-type 00, and the G-type AF phase is accompanied 
by G-type 00»s!i Here we shall consider in more detail 
G-type 00, relevant for the case of LaV03. We will show 
in particular that the superexchange spin-orbital modeffi 
allows one to understand the microscopic reasons behind 
the G-AF phase observed in LaV03, and predicts that 
|«7 C | ~ J a b, as actually observed^ The comparable size 
of FM and AF exchange constants J c and J a b, respec- 
tively, is unexpected when considering the Goodenough- 
Kanamori-Anderson rules, 34 which would suggest that 
\J C \ is by a factor Jh/U smaller. 

In the case of the vanadates the superexchange 
between the S = 1 spins of V 3+ ions in a per- 
ovskite lattice results from virtual charge excitations 
{t\ g )i{t\ g )j ^ (*l 9 )i(*2g)j- These charge excitations in- 
volve the Coulomb interactions in the d 3 configuration 
of a V 2+ ion, parametrized for a pair of t 2g electrons, as 
for the titanates, by the intraorbital Coulomb clement 
U, and by Hund's exchange element Jr [see Eq. (|5.3|) 
and Table H]. 39 The excitation spectrum which leads to 
the superexchange model includes three states^ (i) a 
high-spin state 4 A 2 at energy U — 3Jh, (h) two degen- 
erate low-spin states 2 Ti and 2 E at energy U, and (iii) 
a 2 T 2 low-spin state at energy U + 2J H (Fig. P). We 
parametrize it by two coefficients: r± — 1/(1 — 3?y) and 
7-3 = 1/(1 + 2??) . A general Hamiltonian was already 
given in Ref. Il(:t here we shall analyze it assuming that 
the xy (|c)) orbitals are singly occupied at each V 3+ ion, 
as concluded from experiment22iS4 and from electronic 
structure calculations^ Therefore, the electron densi- 
ties in the remaining two orbitals satisfy at each site i 
the local constraint, 

n ia + n ib = 1. (6.1) 

One finds then that the superexchange TLu{d 2 ) for a bond 
(ij) along the c axis consists of 

H[ c) = -^r 1 (2 + S r S J ){^-f l -r J ), (6.2) 

x(^-^;-^;+5-M> (6.3) 
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and for a bond in the ab plane, 
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(6.7) 



where the operators fi describe orbital pseudospins T = 
1/2 defined (for each direction 7) by the orbital doublet 
{\yz), \xz)} = {\a), \b}}. At each site i there is precisely 
one electron in these two orbitals, and both of them are 
active along the c axis. 

It has already been shown before^ that the quan- 
tum fluctuations play a decisive role in LaVC>3, and the 
observed— anisotropy and temperature dependence of 
the high-spin excitations in optical spectroscopy is re- 
produced when the theory includes them. In the next 
Sec. I VI Bl we will present now the limitations of a sim- 
plified approach which is widely accepted^ and assumes 
that the 00 in LaV03 is quite rigid already at the mag- 
netic transition. Quantum fluctuations lead to important 
corrections which go beyond this picture and modify the 
temperature variation of the optical spectral weights, as 
we discuss in Sec. IVI CI 



B. Spin exchange constants and optical spectral 
weights for alternating orbital order 

After decoupling of spin and orbital variables, the ef- 
fective spin exchange constants J c and J ao in Eq. (|2.11|) 
can be obtained from Eqs. H6.2|) - (|6.7|l by averaging over 
orbital correlations, as derived in Ref. |2(J, 



Jr. = 



Jab — 



1 

4 



Wi - ( r i - vn - Ws) 



- 2r?r 3 (rf r. y ) 



1 - 777*1 - 777-3 



+(n - r\r x - r?r 3 )(- + (r?r/) 



(6.8) 



(6.9) 



They are determined by the orbital correlations, which 
result not only from the superexchange Tiu(d 2 ) but 
also from the interactions with the lattice (from the JT 
term)^& as discussed in detail in Sec. IVI CI Here we shall 
consider first classical states with alternating \a) and \b) 
orbitals in the ab planes: G-type 00 when these orbitals 
alternate also along the c axis, and C-type 00 for re- 
peated either \a) or \b) orbitals along the c axis. These 



classical states are frequently assumed as a consequence 
of a strong JT term and observed GdFe03 distortions^ 
Such classical order would naturally follow from a strong 
JT effect, but it is still controversial whether the JT in- 
teraction is actually that strong in the vanadates. Fortu- 
nately, there are already experimental results which can 
help to resolve this controversy, and we address this issue 
in more detail in Sec. IVI CI 

First we consider the limit of strong JT interaction 
with rigid 00. This implies (n-fj) ~ (rf r?) = — \ along 
each bond for G-type 00, so one finds from Eqs. I|6.8() 
and 1)6.9(1 fixed values of the exchange constants: 



j: 



C-AF 



j: 



C-AF 



= \j[l-*l{ri+r z )]. 



(6.10) 
(6.11) 



The FM interaction J^~ AF increases in lowest order lin- 
early with Hund's exchange 77, and the above values of 
the exchange constants give G-AF spin order. 

Whether or not such a classical 00 is realized in the 
ground state, can be investigated by analyzing the con- 
sequences of the present theory for the distribution of 
spectral weight in the optical spectroscopy. The optical 
spectral weights follow by averaging the individual contri- 
butions to the superexchange, see Eqs. (|2.13(l and l|2.15|l . 
One finds for G-type 00 the optical spectral weights: 
— for polarization along the c axis, 



^Jri{s c 



2, 



4 c) = 0, 



(l-s c ), 



for polarization in the ab plane, 



R (ab) = 



K. 



K. 



(ab) 
2 

(ab) 



lj ri (s ab + 2), 
5 

— J(l-s ab ), 
ijr 3 (l - s ab ). 



(6.12) 

(6.13) 
(6.14) 

(6.15) 
(6.16) 
(6.17) 



As in all other cases, they depend on two model parame- 
ters, J and 77, and on the spin-spin correlations {s c , s ao } 
(the 00 is already fixed). 

Again, as in the case of LaMn03 (Sec. IV), the anal- 
ysis of the optical spectra suggests that the effective pa- 
rameters are somewhat different from the atomic values, 
primarily due to the screening of both Coulomb U and 
Hund's exchange Jh- We use here the parameters de- 
duced from the neutron experiments^ and from the op- 
tical spectra™ — one finds^S J = 40 meV and 77 = 0.13. 
These values imply that J c ~ —4.3 meV and J a b — 6.8 
meV, which lead to a MF estimate of the transition tem- 
perature Tj^P ~ 270K. This value has to be still reduced 
by an empirical factor— (close to 68% for S = 1) to 
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estimate the effect of thermal fluctuations, so one finds 
TV ~ 180 K in reasonable agreement with the experi- 
mental value of 140 Ki2i We note that the above values 
of the exchange constants are in good agreement with the 
neutron experiments^ 

It is instructive to test this classical approach by ana- 
lyzing its predictions for the optical spectral weights. We 
evaluated the spin correlation functions for bonds in the 
ab plane and along the c axis using the Oguchi method, 
and used an order parameter (S z ) at neighboring sites, 
which acts on the considered bond by MF terms and 
vanishes at Tn — 0.4J, as explained in Ref. (2(| The 
correlation s c for FM bonds along the c axis can be ob- 
tained analytically — while s a b for the AF bonds was de- 
termined by a numerical approach described in Appendix 
IO As expected, the spectral weight due to the high- 
spin excitations dominates for c polarization. However, 

(c) 

when the 00 is perfect, the anisotropy between K\ 

and K[ ah) at T = 0, being close to 8:1 (Fig. O, is now 
less pronounced than in the case of joint spin and orbital 
dynamics—I At low temperature the low-spin excitations 
dominate the optical spectral weight for ab polarization. 
One finds that the low-energy spectral weight along the 

(c) 

c axis K\ decreases with increasing temperature. Si- 
multaneously, the low-energy spectral weight in the ab 
planes A^ ab ' increases, and the anisotropy goes down to 
- 5 : 2 at T = 0.85J (i.e., at T ~ 300 K for T N = 140 
K). It is quite remarkable that the present classical ap- 
proach gives qualitatively a very similar distribution of 
the spectral weights and their temperature dependence 
for the FM and AF bonds in LaV0 3 and in LaMn0 3 (cf. 
Figs. 151 and H^. 

A crucial test of the present theory concerns the tem- 
perature dependence of the high-spin spectral weight 

(c) 

along the c axis K\ , which according to experiment 
decreases by about 50% between low temperature and 
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FIG. 12: Kinetic energy terms (in units of J) per bond l|2.13fl . 
as obtained for the C-AF phase of the cubic vanadates: high- 
spin (solid lines) and low-spin (dashed lines) and 
(long-dashed lines), for increasing temperature T/J: (a) 
along the FM c axis (7 = c); (b) within the AF ab plane 
(7 = ab). In part (a) the experimental points from Ref. 19 ll for 
the low-energy spectral weight were reproduced, after proper 
rescaling to match the value of at T — > 0. Parameters: 
77 = 0.13. 



TABLE VI: Exchange constants J 7 and optical spectral 
weights -Kl 7 ' (all in units of J) as obtained for: LaV03 (G- 
AF phase with G-type OO), and for YVO3 (G-AF phase with 
G-type OO) at T = and above the magnetic transition at 
T = 0.85 J (~ 300 K). Parameters: 77 = 0.13 and T N = 140 
K. 



LaVOa YVO3 
G-AF phase G-AF phase 



T(K) 





300 





300 


Jc 




-0.107 




0.224 


Jab 




0.171 




0.171 




1.640 


1.181 


0.0 


0.0 




0.0 


0.280 


0.586 


0.494 




0.0 


0.0 


0.465 


0.392 


K (ao) 


0.219 


0.471 


0.249 


0.471 


K (ab) 


0.916 


0.532 


0.871 


0.532 


K (ab) 


0.436 


0.253 


0.415 


0.253 



T = 300 Ks-S In the present theory based upon frozen OO 
this decrease amounts only to about 27% [Fig. IT2T a)]. 

(c) 

and the maximal possible reduction of K\ reached in 
the limit of T — > 00 is by 33%. This result suggests that 
the frozen OO scenario in LaV03 is excluded by exper- 
iment; further arguments supporting this point of view 
are given in Sec. IVI CI Note also that the actual values 
of k[ c \ shown in Refs. |2(] and [3 had to be reduced to 
match the classical prediction at T = 0. 

A second classical state with C-type OO, as proposed 
for the low temperature G-AF phase of YVO^r^ gives an 
AF exchange constant along the c axis, 

Jf- AF = ~ J(l - r?r 3 ), (6.18) 

while the value of J^ b ~ AF is the same as J^ b ~ AF , see Eq. 
(|6.11|) . so in the present classical approach it does not 
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change when the 00 changes from the G-type to C-type. 
In contrast, depending on the type of 00, the exchange 
constant along the c axis can be either FM or AF, as 
suggested by experiment^ and confirmed by the data in 
Table IVII The actual values estimated with the same 
parameters as for LaV03 are: J c — 9.0 and J a b = 6.8 
meV. Here one finds that J c > J a {> as in experiment, but 
the agreement with the experimental values^ J° x p = 3.1 
and J^ p = 2.6 meV, is much poorer. Hence, we conclude 
that quantum effects beyond the present classical analysis 
such as orbital fluctuations^ orbital Peierls effect)22i2S 
and spin-orbit coupling, 100 play here an important role 
and have to be included in a quantitative theory. 

Furthermore, while the weights in the ab planes are 
given by spin-spin correlations s a f,, as in Eqs. Ij6.15|) - 
(|6.17|) . for the spectral weights in the C-type 00 one 
finds along the c axis 

K[ c) = 0, (6.19) 
Kf = \j(l-s c ), (6.20) 

= ljr 3 (l- Sc ). (6.21) 

We evaluated these spectral weig hts (see Table[VT|> for the 
standard value of 7/ = 0.13 and for the spin-spin corre- 
lation function s c determined within the Oguchi method 
for an AF bond along the c axis (Appendix [UJ, stabi- 
lized by the above value of the AF interaction J c , and 
taking again Tjv = 0.4J. It is along this axis where the 
theory predicts a markedly different behavior of the spec- 
tral weights from that found above for G-type 00, see 
Table IVII In spite of AF bonds in all three directions, 
the spectral weight distribution is again anisotropic - 
the weights are considerably higher for the ab planes due 
to the broken cubic symmetry in orbital space. In ad- 
dition, the spectral weights obtained in ab planes at low 
temperature differ between the two AF phases, as the 
MF terms are larger in G-AF phase and make this case 
somewhat closer to the classical limit of s a {, = — 1. Par- 

(c) 

ticularly the prediction that K\ — 0, following from the 
classical G-type 00, is easy to verify. In fact, the exper- 
imental data are puzzling as one finds instead finite and 
temperature dependent spectral weight also for the low- 
energy regime in the G-AF phase of YVO^^^ pointing 
out once again that the present calculation with frozen 
00 is oversimplified. 

C. Composite spin-orbital dynamics in LaVOs 

Finally, we demonstrate that the scenario of strong JT 
interaction, quenching the orbital dynamics, cannot ap- 
ply to LaV03. We do so by investigating its consequences 
for the orbital transition temperature T and for the tem- 
perature dependence of the optical intensity k[ c ^ of the 
lowest multiplet transition. Consider first the transition 
temperature T associated with the phase transtion into 



the state with 00. We have already seen in Sees. Ill 
and IV that strong orbital-lattice coupling in a perovskite 
structure would in fact necessarily decouple orbital and 
spin degrees of freedom and lead to T Q ^ Tn, contra- 
dicting the experiment^ 

The JT coupling between the JT-active local lat- 
tice modes Qi and the pseudospin (orbital) variables r? 
(which refer to the active {|a), \b)} orbitals along the c 
axis) may be written as follows, 

n 3T = g J2QiTi+lYlQl ( 6 - 22 ) 

i i 

where Qi is the appropriate linear combination of coor- 
dinates of the ligand ions next to site i. The above local 
coupling induces local distortions and an associated en- 
ergy gain of -Ejt = g 2 /8 (the JT energy) per site, and 
moreover generates a cooperative JT effect in the follow- 
ing way. As the oxygens are shared by two neighboring 
vanadium ions in the perovskite structure of LaV03, the 
JT distortions Qi at nearest-neighbor sites are not inde- 
pendent from each other. Hence the electron densities 
in | a) /| b) orbitals at two vanadium ions at sites i and 
j in the ab plane (see inset in Fig. I13|l get coupled to 
each other, basically because they depend on the actual 
displacement 6 of the shared oxygen ion. More specif- 
ically, this displacement generates interactions between 
the orbitals^iSi and one finds, taking care of the orthog- 
onality constraint on the Qi variables, 

H ab (ij)=J JT T?T*, (6.23) 

with the interaction constant being given by Jjt = 
2Ai?jx = -jAg 2 . Here the coefficient A is determined by 
the phonon spectrum, viz. by all branches in which the 
local c oord inates Q ,; participatejiSi and A < 1 (see e.g. 
Refs. Il02l and for how to estimate A for the per- 
ovskite lattice). The JT interaction H6.23[l reflects the 
cooperative nature of the JT problem. It favors orbital 
alternation and thus supports the superexchange orbital 
interaction J^ h rS and the orbital model relevant for the 
G-AF phase of LaV03 is then 

n T = Jl • 3 -\) + v *b E T * T h ( 6 - 24 ) 

(ij)\\c ' (ij)\\ab 

with V a b = Jab + Jjt- The first term follows from the 
spin-orbital model, and JJ = Jr\. 

While the superexchange contribution to V a b is small 
(adopting the values of Ref. 13, J T ah ~ 2 meV, i.e., 
^ab <C </J — 33 meV) , one finds that yet V ab > Jc if one 
assumes A = 1, i.e. if one basically identifies the JT inter- 
action with the JT energy, and further accepts the esti- 
mate for Ejt ~ 27 meV given in Ref. 26j. Then the Ising 
term quenches the JJ-driven orbital dynamics and leads 
to a cooperative transition at T Q which locks the orbitals 
in all three directions. However, one should be aware that 
the total energy decrease produced by the JT distortion 
of the lattice, Ej T , as obtained in an ab initio calculation 
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JjT/Jc 



FIG. 13: (Color online) Orbital transition temperature T a 
for increasing JT interaction Jjt in: orbital-only model [Eq. 
I|6.24[l . dashed line], and spin-orbital model [Eqs. 1)6,24^ and 
116.2511 . solid line], compared with the Neel temperature Tjv 
(dashed-dotted line). Both T and Tm in units of JJ. The 
inset shows intersite orbital correlations due to JT distortions 
in LaV0 3 . 



such as in Ref. |26|, actually equals E' JT = Ejt + -Ej T , 
i.e. it comprises both the local energy gain Ejt, which 
does not contribute to the orderingpiSi as well as the JT 
ordering energy Ej T = z a fc Jtt/8 = A-Ejt- The JT contri- 
bution to V a b is therefore given by Jjt = 2Ai? JT /(l + A), 
and since A is usually appreciably smaller than 1 one 
should expect that Jjt is definitely smaller than -E" JT , so 
that it is more likely that actually V a b ^ JJ- 

To gain more insight in the role of the JT interaction 
we performed exact diagonalization of eight-site chains 
along the c axis, combined with a MF treatment of the or- 
bital interactions in the ab plane, to determine the orbital 
transition temperature T . As expected, it scales with 
Jjt in the regime of large JT interaction Jjt > 0.1 JJ 
(but still JjtJS 2 JJ) as shown in Fig. ^| as was also 
found in Ref. 1261 However, we were surprised to see that 
Motome et alm& have discarded this result considering 
it to be "an artefact of the MF treatment" , and argued 
that the coupling between neighboring sites is determined 
solely by JJ . The latter applies only in the limit of ex- 
tremely anisotropic coupling, V a b <C JJ, which is not rel- 
evant here. In fact, it is a classical resul1*i2i that the JT 
coupling of Eq. (|6.22|) also induces intersite interactions 
between the orbitals as in Eq. (|6.23|) which may actually 
dominate over the superexchange term and determine T 
in the limit of large g. As T Q is 20-80% above Tjv in 
other RVO3 (R= Ce, Pr, Nd, Y, etc.) compounds^ it is 
not plausible that for R=La it is determined by superex- 
change alone. Therefore we argue that as the ab plane 
correlations are of Ising-type, the MF result should be in 
fact a reasonable estimate for T . Thus, the proximity 
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FIG. 14: (Color online) Temperature dependence of: (a) op- 
tical spectral weight of the high-spin transition, , in units 
of J (J is the principal energy scale of the spin-orbital model), 
and (b) orbital correlations {ft • ft+i) along the c axis (bot- 
tom), as obtained for — the spin-orbital model Eqs. 16.26L 
as given in Ref. HH with small Jjt = 0.2JJ (solid lines), 
and — the orbital model Eq. 16.241 with Jjt = 1.6JJ, as 
analyzed in Ref. [2^ (dashed lines) . 



T ~ Tjv (— 0.4 J), implies that the JT interaction is 
small, Jjt ~ 0.1JJ, and the JT splitting of the xz/yz- 
doublet is actually smaller than the superexchange en- 
ergy scale itself, A JT = z a fcJjT/2 = 2J JT ~ 0.2JJ. 

Furthermore, as the spin and orbital exchange interac- 
tions are interrelated ) 20 ' 26 spin disorder should reduce the 
effective orbital exchange JJ. Indeed, we have verified 
that this follows from the full spin-orbital superexchange 
model 16 which contains an extra term, 



1 



{ij)\\c 



1 



-v 



(6.25) 

describing the coupling of spins and orbitals along the 
c axis. However, T Q obtained from the complete spin- 
orbital model, 



Hu{d 2 ) — ri T + Tisr, 



(6.26) 



has about the same value as that found from ri T alone 
(see Fig. and the estimate Ajt — 0.2JJ remains 

valid. 

Moreover, one finds that for large Jjt (e.g., Jjt = 
2E' JT = 54 meV, i.e. for A = 1 and E' JT = 27 meV ~ 



of the orbital and the magnetic transitions in LaVC^^ 



0.8JJ, following Ref. 26) the temperature dependence of 

(c) 

the optical intensity K\ derived from ri T is weak below 
300 K (Fig. I14fl . so adding orbital correlations in this 
range of parameters to the spin correlations of Sec. IVI Bl 
cannot improve the agreement with experiment. Hence, 
this analysis clearly shows that: (i) a substantial vari- 
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ation of K\ c) below 300K (~ J, see Fig. O observed 
in Ref. I29I is obtained only for small Jjt; (ii) while nei- 
ther spin correlations for frozen 00, nor orbital correla- 
tions that follow from 7i T alone would suffice, only the 
full spin-orbital model (|6.26() that includes coupled spin- 
orbital fluctuations is able to explain a large enhancement 

(c) 

of K\ at low temperature. This was indeed demon- 
strated in Ref. I20L where spin and orbital correlations 
were treated self-consistently, and only then the strong 

(c) 

temperature dependence of K{ ' could be successfully re- 
produced by the theory. 

Summarizing, on comparing the results of Figs. 1 131 and 
Hlwith the experimental data)** one has to conclude that 
the proximity of T Q and Tjv in LaV03 implies that Jjt 
is in fact much smaller than the total energy associated 
with the cooperative JT distortion E' JT ~ 27 meV es- 
timated in Ref. l26i The enhanced optical conductivity 
along the c axis2£ also suggests that the local splitting 
Ajt is smaller than the dynamical orbital exchange JJ, 
thus supporting the scenario of fluctuating orbitalsi& in 
LaV03. Therefore, the assumption of rigid 00, which 
was so successful for LaMn03, fails for LaV03 and the 
full quantum spin-orbital many-body problem has to be 
treated explicitly^ 

VII. SUMMARY AND CONCLUSIONS 

The main purpose of this paper was to make the exper- 
imental consequences of the superexchange spin-orbital 
models for correlated transition metal oxides with orbital 
degeneracy more transparent. We formulated a general 
approach to the spectral weights in optical spectroscopy 
and illustrated it on several examples with different mul- 
tiplet structure. While a general feature of all the su- 
perexchange spin-orbital models is a tendency towards 
enhanced quantum fluctuationsj2ii£ we gave reasons why 
in many situations such fluctuations are quenched. One 
then arrives at much simpler reduced models, where cer- 
tain states with 00 allow for a good insight into the 
mechanisms responsible for the magnetic interactions and 
for the optical spectral weights. The common feature of 
all these cases is that the knowledge of only a few ef- 
fective parameters, the superexchange energy J, Hund's 
exchange rj and the charge-transfer parameter R, is suf- 
ficient to work out the quantitative predictions of the 
theory for a given type of orbital ordered state. In some 
of these cases the theoretical models simplify so much 
that it is even possible to perform calculations with the 
help of a pocket calculator. 

The cases of copper fluoride KCUF3 and the mangan- 
ite LaMn03 turned out to be simpler, and could be 
understood with frozen 00 and quenched orbital dy- 
namics below the structural transitions which occur at 
much higher temperatures than the Neel temperature, 
T s 3> Tat. However, we have also seen that particularly 
in tig systems, in the cubic titanates and vanadates, the 
orbital dynamics may not be quenched. Therefore, in 



some cases only the full quantum many-body problem 
gives proper answers for the experimental situation. 

We came to these conclusions by analyzing in detail 
the consequences of decoupling of the spin and orbital 
degrees in states with rigid 00 and by comparing the 
predictions of the theory with the experimental data, 
wherever available. In the undoped manganite LaMn03 
we could provide a consistent explanation of the mag- 
netic and optical experimental data by deducing the val- 
ues of the above effective parameters {J, r\, R}, and next 
showing that both the magnetic exchange constants J a b 
and J c , and the anisotropy and the temperature depen- 
dence of the low-energy optical spectral weights can be 
reproduced by the theory in a satisfactory way. In the 
case of the copper fluoride KCUF3 optical data were not 
available, but the constraints in the theory given by the 
exchange constants are so strong that we could conclude 
that the insulating state in this compound has charge- 
transfer character. It remains to be verified by future 
experiments to what extent the predictions made here 
concerning the optical spectral weights and based on the 
classical picture with ordered e g orbitals apply to KCUF3. 

Also for LaTi03 and YTi03 we investigated the classi- 
cal states with 00 given by certain wavefunctions which 
guarantee that the observed isotropic AF or FM states 
are realized. As in all other cases, the theory predicts in 
such states a rather pronounced temperature dependence 
and spectral weight transfers in the optical spectra near 
the magnetic transition. Future experiments will have 
to establish whether and to what extent such a scenario 
relaying on rigid 00 could be valid. However, already 
without these data we could demonstrate, by looking at 
the exchange constants, that there are certain indications 
that orbital fluctuations play a role and thus the quan- 
tum physics might dominate here over the thermal fluc- 
tuations of the spins alone. 

The case of the cubic vanadate LaV03 really shows 



TABLE VII: Values of the effective parameters of spin-orbital 
models: superexchange J (in meV), Hund's parameter r/, the 
CT parameter 7?, and the microscopic parameters consistent 
with these effective parameters: intraorbital Coulomb inter- 
action U, Hund's exchange Jh, and the energy of the lowest 
CT excitation A (all in eV), deduced from the present anal- 
ysis of the magnetic and optical properties of representative 
transition metal compounds with perovskite structure. The 
values of Jh and U in case of LaVOs and LaMnOs were ob- 
tained from the optical spectra, while the ones for KCUF3 are 
the same as in Ref. I5H 
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that one may encounter the full complexity of the spin- 
orbital superexchange model when the spins and orbitals 
fluctuate coherently and these fluctuations are essential 
to get a meaningful quantitative description of the opti- 
cal data. This case in particular demonstrates the im- 
portance of combining the magnetic and optical data. 
Whereas the classical analysis of the spin fluctuations for 
frozen orbitals seems to suffice to explain the exchange 
constants in LaVC>3, by looking at the optics one realizes 
that the picture of frozen 00 induced by large JT cou- 
pling is here misleading, and the complete spin-orbital 
dynamics has to be considered instead. 

Before concluding this paper, we summarize in Table 
IVIII the effective parameters of the spin-orbital modelsjiS 3 - 
and the possible values of the microscopic parameters — 
the Coulomb interaction U, Hund's exchange Jh, and 
the charge transfer energy A, that are consistent with 
these effective parameters. The values of the superex- 
change constant are 33 < J < 40 meV for KCUF3 and 
LaV03, while for LaMn03 the value of this parameter 
is much higher, J ~ 140 meV. This difference reflects 
a rather high value of the effective hopping parameter t 
in the undoped manganite. By considering the optical 
data of LaV03 and LaMn03 we came to the conclusion 
that Hund's exchange Jh is somewhat reduced from the 
respective atomic values, 54 and we give already these re- 
duced parameters in Table IVTT1 accompanied by the cor- 
responding values of U. 

Summarizing, we have illustrated a common approach 
to the optical and the magnetic data in Mott insula- 
tors with orbital degeneracy, which provides the basis 
for a better theoretical understanding of the experimen- 
tal constraints on the underlying phenomena and on the 
model parameters. It is a unique feature of these sys- 
tems that the superexchange interactions, and the spin, 
orbital, and composite spin-and-orbital correlations in- 
duced by them, are responsible for the distribution of 
spectral weight in the optical excitations. We hope that 
extending the present analysis of the high-energy excita- 
tions in the upper Hubbard band by an analysis of the 
low-energy excitations that occur in doped systems, will 
allow to develop a quantitative theoretical approach de- 
signed to describe the optical spectra of doped transition 
metal insulators. 
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APPENDIX A: EFFECTIVE MODEL FOR 
LaMn0 3 

Here we present an analysis of the magnetic and op- 
tical data within an effective d — d model for LaMn03, 
given by the Hu{d 4 ) term alone, i.e., assuming R = 
in Eq. (|4.1U|) . Unlike in KCuF 3 , the alternating OO in 
LaMn03 refers to the orbitals occupied by electrons, and 
is characterized by a single angle 8 in Eqs. I|3.17|l . We 
performed an analysis of the exchange constants J c (|4.23|) 
and J b 14.22[) as functions of the effective parameters J 
and rj. Although these parameters cannot be uniquely de- 
termined, we have verified that only a narrow range of the 
orbital angle 8 ~ 90° — 100° gives reasonable agreement 
with experiment. Here we present the results obtained 
with J = 170 meV and rj = 0.16. 

Both exchange constants exhibit a rather strong depen- 
dence on the orbital angle 8 (Fig. 115(1 . In contrast to the 
case of KC11F3 , the effective model parameters obtained 
at R = suffices to explain even almost quantitatively 
the observed exchange constants in LaMn03. This re- 
sult is also consistentiSi with the earlier analysis of Ref. 
IT3I Due to the strong dependence of the exchange con- 
stants J a b and J c on the angle 8, one can exclude the OO 
of alternating directional (3a; 2 — r 2 )/(3y 2 — r 2 ) orbitals, 
obtained with 8 = 120° in Eqs. (|3~T7|) . 

While it is frequently assumed that the ti g superex- 
change (|4.20l) is large, the present analysis shows that 
a consistent description of the magnetic properties re- 
quires instead a rather moderate value of Jt in LaMnOs. 
For the present parameters one finds Jt = 1.70 meVjiS 5 - 
and Jt increases by ~ 35% when the CT terms are in- 
cluded (see Sec. II V C|) . Although it might be argued that 
J t ~ 2 meV is too large as the experimental value of T/v 
in CaMn03 is only 110 Kji2& one deals here de facto with 
different values of the U, Jh and A parameters, namely 
with those for Mn 2+ ions instead of the ones for Mn 3+ 
ions, which apply in the intermediate excited states con- 
tributing to the superexchange in CaMn03. Yet the dif- 
ferences in these parameter values cannot be large from 
the very nature of their physical origin, and so Jt for 
LaMn03 cannot differ by more than a factor of 2 from 
the value of the t2 g -md\iced superexchange in CaMnOsji 3 - 
This supports our finding that a small ratio Jt ~ 4x 10~ 3 i 
corresponds to realistic parameter values for LaMnOs . 

Next we consider the temperature dependence of the 
spin-spin correlations and the optical spectral weights 
(|2.15|) . As in the full model discussed in Sec. II V PI 
one may assume frozen OO in the relevant temperature 
range below room temperature, and derive the temper- 
ature dependence of the optical spectral weights from 
that of the intersite spin-spin correlations (|2.12l) . with 
the latter determined in a cluster method as explained in 
Appendix [5] Using the exchange interactions obtained 
with the present parameters at 8 — 94° (Fig. I15fl . one 
finds the empirical estimats^ 3 . T/v — 146 K, which rea- 
sonably agrees with the experimental value T™ p = 136 
K. 62 The AF bonds feel staggered MF terms for T < T N , 
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FIG. 15: (Color online) Superexchange interactions — J a t 
14.221 and J c 14.231 as functions of the orbital angle 9 (solid 
lines), obtained within the effective superexchange model 
H.u(d 4 ) (R = 0). Experimental values^- of —J a b and J c for 
LaMnOa (indicated by circle and diamond) are nearly repro- 
duced for the OO with 9 = 94°. Parameters: J = 170 meV, 
V = 0.16. 
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so the spin-spin correlations s c had to be determined by 
a numerical solution, as explained in Appendix [BJ Both 
correlation functions change fast close to T/v, reflecting 
the temperature dependence of the Brillouin function, 
and remain finite at T ^> Tn [Fig. IPoT a)]. 

The effective model also allows one to discuss the qual- 
itative features of the spectral weight distribution in the 
optical spectra. The theory predicts that only high-spin 
optical excitations are allowed at T = for the FM bonds 
in the ab planes, and one finds for this polarization a 
large kinetic energy [Fig. Iltjf bl]. In contrast, the 

optical excitations for the AF bonds along the c axis are 
predominantly of low- spin character, and thus the ki- 
netic energy is rather small, resulting in a large 

anisotropy K'f^ : K[ c ' ~ 10 : 1 of the low-energy optical 
intensities confirmed by experiment. 28 When the tem- 
perature increases and the spin-spin correlations weaken, 
this anisotropy is reduced, but remains pronounced also 
at T > T N and still exceeds 2 : 1 at T = 300 K due to 
the persisting OO. As shown elsewhere^ also a quantita- 
tive analysis of the present effective model gives a rather 
satisfactory agreement with the optical data— in the en- 
tire temperature range. In contrast, the total optical 
intensities have a much weaker temperature dependence 
and anisotropy [Fig. llbT b)]. Thus, the main features on 
the experimentally observed intensity distribution and its 
temperature variation in the optical spectra— are well re- 
produced already by the present effective model. 

Finally, we verify whether the used parameters J = 170 
meV and r\ = 0.16 could be derived from the microscopic 
parameters of the CT model. The value of Jh = 0.90 



FIG. 16: (Color online) Temperature dependence of: (a) spin- 
spin correlations {Si • 5j) 7 for 7 = a, c, and the order parame- 
ter (S z ) (dashed line); (b) kinetic energies K[ l] for t he hi gh- 
spin excitations, standing for the optical intensities 12.151 at 
low energy along 7 = a,c axes (solid lines), and total kinetic 
energies (long-dashed lines), as obtained for LaMnOa 

within the effective superexchange model H.u(d 4 ) for the OO 
given by an angle 9 ~ 94° . Parameters as in Fig. 1151 



eV follows from the spectroscopic values of the B and C 
Racah parameters^ so rj — 0.16 implies U ~ 5.6 eV. 
Knowing the value of J, this leads to an estimated effec- 
tive d— d hopping element t ~ 0.49 eV. Indeed, one finds 
that these microscopic parameters are in the expected 
range. The effective Coulomb interaction was estimated 
within the effective d — d model for LaMnOa as U ~ 5.5 
eV from spectroscopic data^ so the agreement is close to 
perfect. Furthermore, taking the usually accepted values 
of t p d = 1.5 eV and A = 5.0 eV, following Refs. and 
I43I one finds a very plausible value of the effective hop- 
ping parameter t = t^ d / A = 0.45 eV, again quite close 
to the value derived above. Note that the experimen- 
tal magnetic interactions in doped bilayer manganites 
were explained with a similar value of t = 0.48 eV fixed 
by experiments^ This favorable comparison emphasizes 
once again our main conclusion that the relevant model 
parameters can be derived by combining the results of 
magnetic and optical experiments, whenever available. 
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APPENDIX B: SPIN-SPIN CORRELATIONS IN 
LaMn0 3 

Here we describe briefly a simple method which we 
used to determine the spin-spin correlation functions s ab 
and s c for a pair of interacting spins S — 2 at nearest 
neighbor Mn ions in LaMn03. The spin- spin correlations 
were obtained by performing a statistical average over the 
exact eigenstates for a single (FM or AF) bond, found in 
the presence of MF terms acting at each spin of the bond 
and originating from its neighboring spins. This, in fact, 
is the simplest cluster MF theory, known in the theory 
of magnetism as the Oguchi method^ As long as the 
MF terms vanish (at T > Tjv), one finds for the various 
eigenstates with degeneracy (2Stot + 1), labeled by the 
total spin Stot = 0, 1, • • • , 4, 



x exp(-4/3|J afc |), (B7) 

z 2 = (1 + 2cosha; + 2cosh2x)exp(-7/3|J afc |), (B8) 

zi = (1 + 2cosha;)exp(-9/3|J Q6 |), (B9) 

z = cxp(-10/3|J afc |), (BIO) 



with x = (3h ab and f3 = 1/fceT. Note that the term 
oc Z3 is absent in the numerator of Eq. (|B5|I . because 
(Stot = 3|& • SjlStot = 3) = [see Eq. (EE)]. 

For an AF bond (ij) along the c axis the Hamiltonian 
is given by 



H^f = jJi ■ S* - h c (S? 



where the molecular field 

h c = (4\J ab \+J c )(S z 



(Bll) 



(B12) 



(S'tot I Si ■ Sj I S'tot/ 



r S'tot (St, 



1 



G. 



(Bl) 



Depending on S'tot the scalar product (Stot|S, • Sj|Stot) 
varies between —6 and 4. 

Consider first a FM bond (ij) in the ab plane, with the 
Hamiltonian given by, 

H\f = -\J ab \S t ■ Sj - h ab (S z + S?), (B2) 

where for the A-AF phase with uniform order parameter 
in the ab plane (S z ) = (Sf) the identical MF acting at 
both spins is 



h ab = {3\J ab \+2J c )(S z ). 



At T > Tjv the eigenenergies of follow from Eq. 

HBljl. For T < Tv the order parameter (S z ) could in 
principle be determined self-consistently in the present 
cluster approach. However, to gain a qualitative insight 
into the temperature dependence of s ab it suffices to use 
a self-consistent solution of the MF equation, 



(B3) 



alternates between the sites i and j, so (Sf) = 
-(S*|) = (S z ). Unlike for the FM bond, the field 
h c couples now states which belong to different val- 
ues of Stot (but to the same value of S t 2 ot ). The 
staggered MF plays no role for Sf ot = 4, while for 
S z ot < 4 the eigenstates have been found by di- 
agonalizing the matrices, with diagonal elements fol- 
lowing from Eq. ijBlfl . and offdiagonal ones oc h c : 
- for SfL = 3, 



for S t 2 ot - 2, 



(B13) 




4J C 



^ h 



for Si 



f h c 



-4 



-3,L 



(B14) 



(S z ) = 



2 2T 



(S z 



25+1 , (2S+1T N/CIZ , 

coth ' 

2 

2 coth (--(5 



(B4) 



As the MF value fc B T^ IF = 4(2|J ah | + J c ) is overesti- 
mated, it is appropriate to use in Eq. I|B4(I the value 
of Tjv after an empirical reduction^ leading to Tv ~ 
0.705T^ F . One finds an analytic solution for s ab - 



Sab 



4^4 — 3^2 — 5z\ — 6zq 

Zi + Z3 + Z 2 + Z\ + Zq ' 



(B5) 



where the terms z% refer to the subspaces of total spin 

Stot, 

z 4 = 1 + 2 cosh x -I- 2 cosh 2x + 2 cosh 3a; 

+ 2cosh4x, (B6) 

Z3 = (1 + 2 cosh x + 2 cosh 2x + 2 cosh 3a;) 
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In this way a complete set of eigenstates {\n)} with en- 
ergies {E n } for n = 1, 2, • • • , 25 was determined. Finally, 
the spin-spin correlation function s c was found using a 
standard formula 

Sc = \ ■ ^» exp(-/?£„), (B17) 



where Z = J2 n exp(— j3E n ) is the partition function. 



Here 



h c = (4J ab + \J c \)(S z ), 



(C6) 



x — (3h c and (3 = compare with Eq. 1B12|1 . 

For an AF bond {ij) in the ab plane the MF Hamilto- 
nian is given by 



H 



(ab) 



— JabSi • Sj h a b(Si ^jn 



(C7) 



APPENDIX C: SPIN-SPIN CORRELATIONS IN 
LaV0 3 AND IN YV0 3 

The short-range spin-spin correlations s 7 for the cubic 
vanadates were determined using the Oguchi method^ 
for a bond of interacting 5 — 1 spins. As in the case 
of LaMnC>3 (see Appendix [BJ) , we solve exactly a single 
FM (AF) bond {ij) with interaction J c (J a b), and the 
MF terms oc (5 2 ) originating from neighboring spins and 
acting on each spin of the bond. In the present case the 
scalar product is given by 

{Stat\Si • S,|Stot> = ~5 tot (5 tot + 1) - 2. (CI) 

For a FM bond, now along the c axis, one finds an 
analytic solution^ This problem is analogous to that 
given by Eq. I|B3(1 . Using the MF approximation, the 
order parameter (5 2 ) was determined from Eq. 1|B4[1 with 
5 = 1, and T™ F = A(2J ah + |J c |)/3 was reduced to T N ~ 
0.684T^ F as appropriate for 5 = 1 spins. 73 The final 
result for s r reads 



Z2 - z\ - 2z Q 



Z2 + Z\ 



(C2) 



where the terms Zi originate from different subspaces of 
total spin 5 to t = 2,1,0, 



Z2 = 1 + 2 cosh a; + 2 cosh 2a;, 
Z\ = (1 + 2 coshx) exp(— 2/3| J c 
z = exp(-3/3| J c |). 



(C3) 
(C4) 
(C5) 



In analogy to an AF bond in LaMn03 (see Appendix lB|l . 
the correlation function s a b can be found numerically by 
considering the subspaces of 5^ ot . The molecular field, 



(3J ab + 2\J c \){S z ) 



(C8) 



does not contribute to the highest eigenenergies Eg t g = 
Jab in the subspace of 5 t z ot = 2 [cf . with Eq. l|B3j) ] , while 
the remaining eigenstates had to be found by diagonal- 
izing the matrices corresponding to other values of Sf ot : 
-for 5^ = 1, 



for = 0, 



(C9) 
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By solving the above eigenproblems, we determined a 
complete set of eigenstates {|n)}, with energies E n , la- 
beled by n = 1, 2, • • • , 9. Therefore, the spin-spin corre- 
lation function s a b for two interactings 5=1 spins on 
an AF bond follows in the present case from an equation 
similar to Eq. I|B17I) . with the relevant matrix elements 
(n|5j • Sj\n) now given by Eq. (jClfl . 



1 M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

2 E. Dagotto, T. Hotta, and A. Moreo, Phys. Rep. 344, 1 
(2001); E. Dagotto, New J. Phys. 7, 67 (2005). 

,j E. Dagotto, Nanoscale Phase Separation and Colossal 
Magnetoresistance, Springer Series in Solid State Sciences 
Vol. 136 (Springer- Verlag, Heidelberg, 2003). 

4 S. Maekawa, T. Tohyama, S. E. Barnes, S. Ishihara, W. 
Koshibae, and G. Khaliullin, Physics of Transition Metal 
Oxides, Springer Series in Solid State Sciences Vol. 144 
(Springer- Verlag, Heidelberg, 2004). 

5 Y. Tokura and N. Nagaosa, Science 288, 462 (2000); Y. 
Tokura, Physics Today, July 2003, p. 50. 



6 P. W. Anderson, Phys. Rev. 115, 2 (1959). 

7 K. I. Kugel and D. I. Khomskii, Usp. Fiz. Nauk 136, 621 
(1982) [Sov. Phys. Usp. 25, 231 (1982)]. 

8 M. Cyrot and C. Lyon-Caen, J. Phys. (Paris) 36, 253 
(1975); S. Inagaki, J. Phys. Soc. Jpn. 39, 596 (1975); C. 
Castellani, C. R. Natoli, and J. Ranninger, Phys. Rev. B 
18, 4945 (1978); 18, 4967 (1978); 18, 5001 (1978). 

9 L. F. Feiner, A. M. Oles, and J. Zaanen, Phys. Rev. Lett. 
78, 2799 (1997). 

G. Khaliullin and V. Oudovenko, Phys. Rev. B 56, 
14243(R) (1997); L. F. Feiner, A. M. Oles, and J. Zaanen, 
J. Phys.: Condens. Matter 10, L555 (1998); G. Khaliullin 
and R. Kilian, ibid. 11, 9757 (1999). 



33 



11 A. M. Oles, L. F. Feiner, and J. Zaanen, Phys. Rev. B 61, 
6257 (2000). 

12 S. Ishihara, J. Inoue, and S. Maekawa, Phys. Rev. B 55, 
8280 (1997); R. Shiina, T. Nishitani, and H. Shiba, J. 
Phys. Soc. Jpn. 66, 3159 (1997). 

13 L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 

14 A. Weisse and H. Fehske, New J. Phys. 6, 158 (2004). 

15 G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 3950 
(2000); G. Khaliullin, Phys. Rev. B 64, 212405 (2001). 

16 G. Khaliullin, P. Horsch, and A. M. Oles, Phys. Rev. Lett. 
86, 3879 (2001). 

17 H. F. Pen, J. van den Brink, D. I. Khomskii, and G. A. 
Sawatzky, Phys. Rev. Lett. 78, 1323 (1997); R. Shiina, 
F. Mila, F.-C. Zhang, and T. M. Rice, Phys. Rev. B 63, 
144422 (2001); S. Di Matteo, N. B. Perkins, and C. R. Na- 
toli, ibid. 65, 054413 (2002); S. Di Matteo, G. Jackeli, C. 
Lacroix, and N. B. Perkins, Phys. Rev. Lett. 93, 077208 

(2004) ; S. Di Matteo, G. Jackeli, and N. B. Perkins, Phys. 
Rev. B 72, 024431 (2005); Phys. Rev. B 72, 020408(R) 

(2005) . 

18 A. M. Oles, Acta Phys. Polon. B 32, 3303 (2001); Phys. 
Status Solidi (b) 236, 281 (2003); ibid. 242, 963(E) 
(2005). 

19 P. Fazekas, Lecture Notes on Electron Correlation and 
Magnetism (World Scientific, Singapore, 1999). 

20 G. Khaliullin, P. Horsch, and A. M. Oles, Phys. Rev. B 
70, 195103 (2004). 

21 J. S. Lee, M. W. Kim, and T. W. Noh, New J. Phys. 7, 
147 (2005). 

22 D. Baeriswyl, J. Carmelo, and A. Luther, Phys. Rev. B 
33, 7247 (1986). 

23 P. Horsch and W. Stephan, Phys. Rev. B 48, 10595 
(1993); H. Eskes and A. M. Oles, Phys. Rev. Lett. 73, 
1279 (1994); H. Eskes, A. M. Oles, M. Meinders, and W. 
Stephan, Phys. Rev. B 50, 17980 (1994). 

24 M. Aichhorn, P. Horsch, W. von der Linden, and M. 
Cuoco, Phys. Rev. B 65, 201101(R) (2002). 

25 J. Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 

26 Y. Motome, H. Seo, Z. Fang, and N. Nagaosa, Phys. Rev. 
Lett. 90, 146602 (2003). 

27 K. Tobe, T. Kimura, Y. Ohimoto, and Y. Tokura, Phys. 
Rev. B 64, 184421 (2001). 

28 N. N. Kovaleva, A. V. Boris, C. Bernhard, A. Kulakov, A. 
Pimenov, A. M. Balbashov, G. Khaliullin, and B. Keimer, 
Phys. Rev. Lett. 93, 147204 (2004). 

29 S. Miyasaka, Y. Okimoto, and Y. Tokura, J. Phys. Soc. 
Jpn. 71, 2086 (2002). 

30 A. A. Tsvetkov, F. P. Mena, P. H. M. van Loosdrecht, D. 
van der Marel, Y. Ren, A. A. Nugroho, A. A. Menovsky, I. 
S. Elfimov, and G. A. Sawatzky, Phys. Rev. B 69, 075110 
(2004). 

31 J. S. Lee, Y. S. Lee, T. W. Noh, S. J. Oh, J. Yu, S. 
Nakatsui, H. Fukazawa, and Y. Maeno, Phys. Rev. Lett. 
89, 257402 (2002). 

32 J. Zaanen, G. A. Sawatzky, and J. W. Allen, Phys. Rev. 
Lett. 55, 418 (1985). 

33 We note, however, that frequently the average energy of a 
CT type excitation A is considered instead and compared 
with U . 

34 J. B. Goodenough, Magnetism and the Chemical Bond 
(Interscience, New York, 1963). 

35 J. Zaanen and A. M. Oles, Phys. Rev. B 37, 9423 (1988); 
H. Eskes and J. H. Jefferson, ibid. 48, 9788 (1993). 

36 M. V. Mostovoy and D. I. Khomskii, Phys. Rev. Lett. 92, 



167201 (2004). 

37 O. K. Andersen, N. Klose, and H. Nohl, Phys. Rev. B 17, 
1209 (1978). 

38 J. Zaanen and A. M. Oles, Phys. Rev. B 48, 7197 (1993). 

39 J. S. Griffith, The Theory of Transition Metal Ions (Cam- 
bridge University Press, Cambridge, 1971). 

40 A. M. Oles, Phys. Rev. B 28, 327 (1983). 

41 Fractional contributions due to exchange terms tx Jh 
shown in the spectrum for the d 5 case [see Fig. 0[a)] fol- 
low from the differences between the exchange elements 
J a (3 in Hamiltonian Hi nt 12. 211 (see Table I), and were ob- 
tained using the relation C ~ 4_B satisfied approximately 
by the experimental values for Mn 2+ ionsi— 

42 D. van der Marel and G. A. Sawatzky, Phys. Rev. B 37, 
10674 (1988). 

43 T. Mizokawa and A. Fujimori, Phys. Rev. B 54, 5368 
(1996). 

44 M. V. Mostovoy and D. I. Khomskii, Phys. Rev. Lett. 89, 
227203 (2002); A. J. W. Reitsma, L. F. Feiner, and A. M. 
Oles, New J. Phys. 7, 121 (2005). 

45 P. F. Maldague, Phys. Rev. B 16, 2437 (1977). 

46 A. Chattopadhyay, A. J. Millis, and S. Das Sarma, Phys. 
Rev. B 61, 10738 (2000). 

47 K. H. Ahn and A. J. Millis, Phys. Rev. B 61, 13545 (2000) 

48 F. Mack and P. Horsch, Phys. Rev. Lett. 82, 3160 (1999) 

49 S. Kadota, I. Yamada, S. Yoneyama, and K. Hirakawa, J 
Phys. Soc. Jpn. 23, 751 (1967). 

50 S. K. Satija, J. D. Axe, G. Shirane, H. Yoshizawa, and K 
Hirokawa, Phys. Rev. B 21, 2001 (1980). 

51 L. Paolasini, R. Caciuffo, A. Sollier, P. Ghigna, and M 
Altarelli, Phys. Rev. Lett. 88, 106403 (2002); R. Caciuffo 
L. Paolasini, A. Sollier, P. Ghigna, E. Pavarini, J. van den 
Brink, and M. Altarelli, Phys. Rev. B 65, 174425 (2002) 

52 M. T. Hutchings, E. J. Samuelsen, G. Shirane, and K 
Hirakawa, Phys. Rev. 188, 919 (1969). 

53 De facto, the values of r^ and r% differ and depend on 
bond direction below the structural transition when the 
JT distortions remove the degeneracy of the e g orbitalsi- 
Here the idealized case of r^ = rz is discussed. 

54 J. Zaanen and G. A. Sawatzky, J. Solid State Chem. 88, 
8 (1990). 

55 A. I. Lichtenstein, V. I. Anisimov, and J. Zaanen, Phys. 
Rev. B 52, 5467(R) (1995). 

56 J. E. Medvedeva, M. A. Korotin, V. I. Anisimov, and A. 
J. Freeman, Phys. Rev. B 65, 172413 (2002). 

57 D. A. Tennant, T. G. Perring, R. A. Cowley, and S. E. 
Nagler, Phys. Rev. Lett. 70, 4003 (1993); B. Lake, D. A. 
Tennant, and S. E. Nagler, ibid. 85, 832 (2000); D. A. 
Tennant, R. A. Cowley, S. E. Nagler, and A. M. Tsvelik, 
Phys. Rev. B 52, 13368 (1995). 

58 A. Millis, Phys. Rev. B 53, 8434 (1996). 

59 J. B. Grant and A. K. McMahan, Phys. Rev. B 46, 8440 
(1992). 

60 D. C. Mattis, The Theory of Magnetism II (Springer- 
Verlag, New York, 1985). 

61 H. Schulz, Phys. Rev. Lett. 77, 2790 (1996). 

62 K. Hirota, N. Kaneko, A. Nishizawa, and Y. Endoh, J. 
Phys. Soc. Jpn. 65, 3736 (1996); F. Moussa, M. Hennion, 
J. Rodriguez-Carvajal, H. Moudden, L. Pinsard, and A. 
Revcolevschi, Phys. Rev. B 54, 15149 (1996); G. Biot- 
teau, M. Hennion, F. Moussa, J. Rodriguez-Carvajal, L. 
Pinsard, A. Revcolevschi, Y. M. Mukovskii, and D. Shuly- 
atev, ibid. 64, 104421 (2001). 



34 



63 Note that a different definition Jh = 2B + C was used in 
Ref.^|to parametrize the spectrum of Mn 2+ ions, ieading 
to a seemingiy smailer vaiue of Jh- 

64 A. E. Bocquet, T. Mizokawa, T. Saitoh, H. Namatame, 
and A. Fujimori, Phys. Rev. B 46, 3771 (1992). 

65 H. Meskine, H. Konig, and S. Sapathy, Phys. Rev. B 64, 
094433 (2001). 

66 Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, I. Koyama, 
M. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
and Y. Endoh, Phys. Rev. Lett. 81, 582 (1998). 

67 J. van der Brink, P. Horsch, F. Mack, and A. M. Oles, 
Phys. Rev. B 59, 6795 (1999); J. van der Brink, New J. 
Phys. 6, 201 (2004). 

68 Although the OO in LaMnOa is in principle un- 
known, we remark that the C-type OO with symmet- 
ric/antisymmetric combinations of e g orbitals, (|z) ± 
\x})/^/2, is favored by a combination of the orbital su- 
perexchange and orbital interactions induced by the JT 
effect at T >T N! H 

69 J. Rodriguez-Carvajal, M. Hennion, F. Moussa, A. H. 
Moudden, L. Pinsard, and A. Revcolevschi, Phys. Rev. 
B 57, 3189(R) (1998). 

70 We have verified that the amplitude of pairs of two or- 
thogonal orbitals is larger for bonds in the ab planes than 
along the c axis, and thus the FM e 9 term is here stronger 
by a factor close to three. 

71 J. van den Brink and D. I. Khomskii, Phys. Rev. Lett. 
82, 1016 (1999). 

72 T. Oguchi, Progr. Theor. Phys. 13, 148 (1955). 

73 G. S. Rushbrooke and P. J. Wood, Mol. Phys. 1, 257 
(1958); see also: M. Fleck, M. G. Zacher, A. I. Lichten- 
stein, W. Hanke, and A. M. Oles, Eur. Phys. J. B 37, 439 
(2004). 

74 A. M. Oles, P. Horsch, and G. Khaliullin, Phys. Status 
Solidi (b) 242, 384 (2005). 

75 M. Raczkowski and A. M. Oles, Phys. Rev. B 66, 094431 

(2002) . 

76 M. Mochizuki and M. Imada, Phys. Rev. Lett. 91, 167203 

(2003) ; New. J. Phys. 6, 154 (2004). 

77 M. Cwik, T. Lorenz, J. Baier, R. Miiller, G. Andre, F. 
Bouree, F. Lichtenberg, A. Freimuth, R. Schmitz, E. 
Muller-Hartmann, and M. Braden, Phys. Rev. B 68, 
060401 (2003). 

78 I. V. Solovyev, Phys. Rev. B 69, 134403 (2004). 

' 9 E. Pavarini, Y. Yamasaki, J. Nuss, and O. K. Andersen, 
New J. Phys. 7, 188 (2005). 

80 E. Pavarini, S. Biermann, A. Poteryaev, A. I. Lichten- 
stein, A. Georges, and O. K. Andersen, Phys. Rev. Lett. 
92, 134403 (2004). 

81 G. Khaliullin and S. Okamoto, Phys. Rev. Lett. 89, 
167201 (2002); Phys. Rev. B 68, 205109 (2003). 

82 S. Ishihara, T. Hatakeyama, and S. Maekawa, Phys. Rev. 
B 65, 064442 (2002). 

83 H. Sawada and K. Terakura, Phys. Rev. B 58, 6831 
(1998). 

84 Equations 15.141 and 15. 1511 were simplified using the iden- 
tity ^{r2+2rz) = T2T3 satisfied by the excitation spectrum 
of Ti 2+ ions. 

85 B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M. v. Zim- 
mermann, J. P. Hill, D. Gibbs, Y. Taguchi, and Y. Tokura, 
Phys. Rev. Lett. 85, 3946 (2000). 



se T. Arima, Y. Tokura, and J. B. Torrance,, Phys. Rev. B 
48, 17006 (1993); Y. Okimoto, T. Katsufuji, Y. Okada, 
T. Arima, and Y. Tokura, ibid. 51, 9581 (1995). 

87 The nearly linear behavior of J c is remarkable and one 
can prove it by performing an expansion in S. Note also 
that for 77 = isotropic exchange interactions J a t = Jc 
are recovered at 8 — 1. 

88 M. Itoh, M. Tsuchiya, H. Tanaka, and K. Motoya, J. Phys. 
Soc. Jpn. 68 2 783 (1999). 

89 Equation 15.26H demonstrates that the FM interaction in 
lowest order increases linearly with Hund's exchange r\ — 



it follows from the identity \(r\ 



r 2 



2r)r 1 r 2 - 



90 C. Ulrich, G. Khaliullin, S. Okamoto, M. Reehuis, A. 
Ivanov, H. He, Y. Taguchi, Y. Tokura, and B. Keimer, 
Phys. Rev. Lett. 89, 167202 (2002). 

91 S. Miyasaka, Y. Okimoto, M. Iwama, and Y. Tokura, 
Phys. Rev. B 68, 100406(R) (2003). 

92 B. Keimer, private communication. 

93 A. V. Mahajan, D. C. Johnston, D. R. Torgeson, and 

F. Borsa, Phys. Rev. B 46, 10966 (1992); H. Kawano, 
H. Yoshizawa, and Y. Ueda, J. Phys. Soc. Jpn. 63, 2857 
(1994); Y. Ren, T. T. Palstra, D. I. Khomskii, A. A. Nu- 
groho, A. A. Menovsky, and G. A. Sawatzky, Phys. Rev. 
B 62, 6577 (2000); M. Noguchi, A. Nakazawa, S. Oka, 
T. Arima, Y. Wakabayashi, H. Nakao, and Y. Murakami, 
ibid. 62, 9271 (R) (2000). 

94 G. R. Blake, T. T. M. Palstra, Y. Ren, A. A. Nugroho, 
and A. A. Menovsky, Phys. Rev. Lett. 87, 245501 (2001); 
Phys. Rev. B 65, 174112 (2002). 

95 H. Sawada, N. Hamada, K. Terakura, and T. Asada, Phys. 
Rev. B 53, 12742 (1996). 

96 S. Miyasaka, S. Onoda, Y. Okimoto, J. Fujioka, M. 
Iwama, N. Nagaosa, and Y. Tokura, Phys. Rev. Lett. 94, 
076405 (2005). 

97 S. Miyasaka, T. Okuda, and Y. Tokura, Phys. Rev. Lett. 
85, 5388 (2000). 

98 C. Ulrich, G. Khaliullin, J. Sirker, M. Reehuis, M. Ohl, 
S. Miyasaka, Y. Tokura, and B. Keimer, Phys. Rev. Lett. 
91, 257202 (2003). 

99 J. Sirker and G. Khaliullin, Phys. Rev. B 67, 100408(R) 
(2003); S. Miyashita, A. Kawaguchi, N. Kawakami, and 

G. Khaliullin, ibid. 69, 104425 (2004). 

100 P. Horsch, G. Khaliullin, and A. M. Oles, Phys. Rev. Lett. 
91, 257203 (2003). 

101 G. A. Gehring and K. A. Gehring, Rep. Prog. Phys. 38, 
1 (1975). 

102 B. Halperin and R. Englman, Phys. Rev. B 3, 1698 (1971). 

103 Having only the magnetic data 85 for LaTiOa, we could 
not obtain meaningful values of the effective parameters 
J and 77 in this case, so this compound was not included 
in Table Em 

104 rp^g an gi e q i s here different from that given in Ref. ^| 
as the value of Jt derived there was underestimated by a 
factor of two. 

Note that this value refers to total spin S = 2 of Mn 3+ 
ions and not to t 2g spin S = 3/2, as frequently defined. 

106 E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 

107 A. M. Oles and L. F. Feiner, Phys. Rev. B 67, 092407 
(2003). 



